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PREFACE 


This compendium summarizes the results of the ERTS-1* investi- 
gations conducted by interdisciplinary teams of the Earth Obser- 
vations Division, NASA Lyndon B. Johnson Space Center at Houston, 
Texas. It presents the important features of each investigative 
area# a summary of the work performed, the processing and analysis 
techniques applied# and the significant findings and conclusions 
resulting from the investigations. 

Seven final reports were prepared for submission to the NASA 
Goddard Space Flight Center in compliance with requirements as 
outlined in the task definition for the ERTS-1 Investigation 
(ER-600) . This project was funded and approved for implementation 
by NASA Headquarters in July 1972. 

The efforts reported in this compendium were performed by 
ERTS-1 analysis team members drawn from the Earth Observations 
Division branches and the ERTS Project Office of the support 
contractor# Lockheed Electronics Company# Inc. The investigations 
are documented in the following individual reports. 
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THE ERTS-1 INVESTIGATION (ER-600) : 

A COMPENDIUM OF ANALYSIS RESULTS OF THE UTILITY OF ERTS-1 
DATA FOR LAND RESOURCES MANAGEMENT 

By R. Bryan Erb 
Lyndon B. Johnson Space Center 

1 .0 THE ERTS-1 INVESTIGATION 
1 . 1 SUMMARY 

This report summarizes the results of the ERTS-1 investiga- 
tions conducted by Earth Observations Division personnel at the 
Lyndon B. Johnson Space Center. Conventional image-interpretation 
and computer-aided classification procedures were the two basic 
techniques used in analyzing the data for detecting, identifying, 
locating, and measuring the surface features related to Earth 
resources. Data from the ERTS multispectral scanner were found 
useful for all applications studied. 

1. The Agricultural Team was able to classify crop types 
(i.e., small grains, truck farm crops, grasses, summer fallow) 

to a 95-percent accuracy on large (12 square hectometers (30 acres) 
or more), well-defined fields. Further breakdowns as to crop 
species (wheat, barley, soybeans, oats, corn) reduced the accuracy 
to 70 to 80 percent for single-date observations. 

2. The Coastal Team was able to map salt marsh areas to an 
accuracy range of 89 to 99 percent and barrier island features uo 
an accuracy range of 87 to 100 percent. 

3. The Forest Team, using conventional and computer-aided 
techniques, was able to differentiate as many as 14 classes of 
forest features to accuracies ranging between 55 and 84 percent. 
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4. The Range Team achieved computer classification accura- 
cies of 70 to 80 percent for separating woodland test fields from 
r.onwoodland test fields. 

5. The Land Use/Urban Team found that by processing data 
from only the urbanized area, it was possible to increase classi- 
fication accuracies of certain urban land use categories. 

6. The Signature Extension Team was able to do temporal sig- 
nature extension on large, relatively homogeneous features. 

7. The computer classification programs were generally more 
accurate than conventional image-interpretation techniques for 
homogeneous features that could be differentiated by spectral 
contrasts. 


1 . 2 INTRODUCTION 

In 1971, the Earth Observations Division (EOD) of the NASA- 
Lyndon B. Johnson Space Center (JSC) submitted a proposal to NASA 
Headquarters for using Earth Resources Technology Satellite 1 
(ERTS-1) data for applications in the Houston Area Test Site 
(HATS). In April 1972, a detailed definition of the functional 
and technical aspects of the project was undertaken. The proposed 
investigations (ER-600) were accepted and funded by NASA Headquar- 
ters in July 1972. R. Bryan Erb, Manager of the EOD Applications 
Office, was named Principal Investigator. 

The main effort of the project involved an assessment of the 
usefulness of the ERTS— 1 data for detecting, identifying, locat- 
ing, and measuring features of applications interest. Its scope 
provided the following: 

1. An applications evaluation in major disciplinary areas 
(forestry, range and grazing management, agriculture, coastal and 
estuarine hydrology, and urban land use management) was initiated. 
It addressed features of prime applications interest because of 
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their geographic and economic or cultural importance. These dis- 
ciplines comprised the major elements of a total information 
hierarchy. This investigation also involved a signature extension 
effort, which assessed techniques for making spatial and temporal 
extensions of established classification. 

2. An Agricultural Stabilisation and Conservation Service 
(ASCS) data-utilization evaluation assessed crops in parts of six 
counties located in different sections of the United States at 
three different periods during the growing season to establish the 
utility of temporal analysis of ERTS-1 data. 

3. A HATS land use study included the analysis of a complete 
ERTS-1 scene by conventional image interpretation and a quarter 
scene by both conventional image interpretation and computer-aided 
analysis. The results were compared with baseline data (high- 
altitude aerial photography) to establish a statistical assessment 
of the ERTS capability for land use classification. 

A report of this scope obviously required the assistance and 
resources of many organizations and individuals, whose help the 
author wishes to acknowledge gratefully. A complete list of the 
numerous contributors to this publication is provided in 
appendix A. 

As an aid to the reader, where necessary the original units 
of measure have been converted to the equivalent value in the 
Systeme International d'Unites (SI). The SI units are written 
first, and the original units are written parenthetically 
thereafter. 


1 . 3 OBJECTIVES 


The general objectives of the JSC ERTS-1 Investigation 
(ER-600) were as follows. 

1. To assess the utility of the ERTS-1 remote-sensor data 
for Earth resources application 

2. To determine the individual and combined roles of the 
ERTS-1 spacecraft, aircraft, and ground data-acquisition systems 
in an integrated Earth resources survey program 

These objectives were based on the concept that, to provide 
resource managers with useful applications information, an entire 
system approach must be considered. This approach would include 
sensor platforms, sensors, data-acquisition subsystems, processing 

and analysis subsystems, and the conversion of data into usable 
products. 

The detailed objectives of the investigations were to 
determine 

1* accu racy and precision to which processed ERTS-1 data 

could be used to identify study features, locate their boundaries, 
and measure their areal size 

2. The nature and uniqueness of the spectral characteristics 
of the study features using the four spectral bands of the ERTS-1 
multispectral scanner (MSS) sensing system with and without atmos- 
pheric corrections 

3. The degree to which temporally composed ERTS-1 data 
derived from registering data sets of the same scene at different 
times could be used to aid in the detection, identification, and 
location of study features 
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1.4 ANALYTICAL APPROACH 
1.4.1 Land Use Hierarchy 

The land use classification system proposed in U.S. Geological 

1 . .... 

Survey (USGS) Circular 671 provided the initial hierarchy of 
study features addressed by each of the ERTS-1 analysis teams 
(table 1-1). The USGS hierarchy is a system designed to be com- 
patible with other classification systems as well as with remote- 
sensor data. The system constitutes first and second levels of 
generalized land use categories that are compatible with remote- 
sensor data acquired from both satellites and high-altitude 
aircraft. 

The USGS land use hierarchy was based on the following 
assumptions . 

1. That Level I categories could be differentiated using 
"satellite imagery with very little supplemental information" 

2. That Level II categories could be derived from "high- 
altitude-aircraft and satellite imagery combined with topographic 
maps " 


The Level I categories would be compatible with imagery 
scales of 1:1,000,000 to 1:250,000, and the Level II categories 
would require imagery scales of approximately 1:100,000 to 
1:24,000. The authors of USGS Circular 671 also noted that sup- 
plementary information might be required to resolve certain 
Level II categories in "especially difficult areas." 


^Although USGS Circular 671 is in the process of revision, 
slight changes in the categories or definitions will not alter the 
significant conclusions of this investigation. 


Level I 


Level II 


Urban and built-up land Residential 

Commercial and services 

Industrial 

Extractive 

Transportation, communications, 
and utilities 

Institutional 

Strip and clustered settlement 
Mixed 

Open and other 

Agricultural land Cropland and pasture 

Orchards, groves, bush fruits, 
vineyards, and horticultural 
area 

feeding operations 

Other 

Rangeland Grass 

Savannas (palmetto prairies) 
Chaparral 

Desert shrub 

Forest land Deciduous 

Evergreen (coniferous and other) 
Mixed 

Nonforested wetland Vegetated 

Bare 

Streams and waterways 
Lakes 

Reservoirs 
Bays and estuaries 
Other 

[salt flats ~ 

Beaches 

Land other than beaches 
Bare exposed rock 
Other 

Tundra 

Permanent snow and icefields 



Permanent snow and icefields 
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The USGS proposed land use hierarchy was designed for remote 
sensors in general, not for a specific sensor, and was structured 
for visual imagery interpretation. Therefore, it was anticipated 
that the system would have to be modified to accomplish best uti- 
lization of automatic data analysis. 

For each of the ERTS-1 investigations, conventional image 
interpretation and computer-aided data analysis were used. As a 
result, some modifications were made to the USGS hierarchy to 
facilitate a specific classification approach. In addition, some 
changes were desirable because of ERTS-1 data characteristics or 
local ground-cover conditions in the areas being classified. How- 
ever, most of the analysis teams addressed study features in 
greater detail than Level II of the USGS hierarchy. As a result, 
some hierarchies at Level III and a few at Level IV were 
established to accommodate the detailed information retrievable 
from the ERTS-1 MSS data. The modifications or extensions of the 
USGS hierarchy are discussed in the individual sections that 
follow. 


1 * 4 *2 Data Processing and Interpretation 

Each EOD ERTS-1 analysis team used data derived from the 
ERTS-1 MSS in its investigation. Several options in conventional 
image interpretation, data enhancement, and computer-aided data 
processing were used, modified, or developed at JSC. The follow- 
ing paragraphs describe the form of the MSS data, the instrumenta- 
available for conventional analysis, the computer programs 
and data flow used for data reformatting and analysis, the tempo- 
ral analysis of data sets, and the atmospheric correction of data. 

Standard products .- The ERTS-1 MSS data are available as 
standard products from the NASA Goddard Space Flight Center (GSFC) 
in film imagery and digital data tape formats. A typical ERTS-1 
scene covers an area of 185 by 185 kilometers (100 nautical miles 
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square) . Film is in either 240- or 70-millimeter format, black- 
and-white or color composites, in the form of transparencies or 
prints. Digital data are stored on computer-compatible tapes 
(CCT's). Four standard 1.3-centimeter (0.5-inch) magnetic tapes 
cover one ERTS-1 scene, with 2340 scan lines and 810 elements 
per line on each computer-compatible tape (CCT) . Both data forms 
were used and analyzed by the ERTS— 1 analysis teams using con- 
ventional image-interpretation techniques and computer-aided 
processing technique s . 

Conventional image-interpretation film data processing . - In 
the conventional image-interpretation film data processing, multi 
band ERTS-1 film imagery was displayed on optical and electro- 
optical color-film viewers to produce color enhancements, which 
accentuate features of interest in contrast to the background 
scene. Photointerpreters delineated, interpreted, and plani- 
metered the enhanced features. Ground- truth information derived 
from low-altitude-aircraft photographs was used to evaluate the 
interpretations and areal measurements obtained. 

The two systems used to create color enhancements were the 
2000 multichannel film viewer (MCFV) and an additive color viewer 
printer (ACVP) . Film density slicing and additive color viewing 
were performed on this equipment. Photograph printing devices 
were available on both systems. 

Computer-aided digital data processing .- The two methods of 
computer-aided data processing used to analyze the digital tape 
data were supervised maximum-likelihood classification and non- 
supervised classification using clustering techniques. The first 
method of processing requires information from training fields 
that are defined by the analyst. Then, unknown data are classi- 
fied into the various classes (i.e. , into the classes that have 
been fed into the computer) . When the classification results are 
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displayed, a threshold is specified so that data points that are 
far from the mean of the class are left unclassified. In the 
second method of processing, the unknown data are organized into 
spectrally homogeneous groups (clusters) and classif ication-type 
clustering maps are produced in which the clusters require identi- 
fication and interpretation in a postprocessing analysis. 

Three main computer systems were available at JSC for digital 
data processing: the Uni vac 1108 computers; the Earth Resources 

Interactive Processing System (ERIPS) , which uses IBM 36C/75 com- 
puters; and the Purdue terminal, a remote terminal connected to 
the IBM 360/67 computing facility at Purdue University. Super- 
vised pattern recognition software packages existed on these three 
facilities. Two clustering algorithms were used for nonsupervised 
classification, the ISOCLS program of JSC and the NSCLAS program 
of Purdue University. (Appendix B contains amplifying information 
about these programs.) Figure 1-1 is a schematic drawing of the 
data flow. 


1.4.3 Other Specialized Analysis Activities 

Additional activities of the analysis teams are described in 
the following paragraphs. 

Conventional image interpretation of first-generation co lor 
composites from digital data .- Two data analysis stations (DAS's) 
at JSC were used to screen and edit digital data. Instead of 
using film imagery to create color enhancements on film viewers, 
the film recorder on the DAS was used to create first— generation 
color composites that were interpreted and analyzed. The positive 
transparency color composites obtained directly from the digital 
data at JSC were of higher quality than those obtained from the 
NASA data-processing facility at GSFC. However, a skewing dis- 
tortion inherent in the digital data display required special 
treatment to restore geometric fidelity. 
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Figure 1-1.- The JSC ERTS MSS data-processing flow. Computer terms, abbreviations, 

and acronyms are defined in appendix B. 
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T emporal analysis .- The cyclic coverage of ERTS-1 made pos- 
sible a temporal analysis of the data. This analysis was per- 
formed to investigate the possible advantages in data 
classification and the analysis gained by a knowledge of the 
changes that occurred, for example, between different vegetative 
phenological stages. The term temporal analysis in this report 
connotes the analysis of a composed data set, which is a unique 
set of data derived from considering a feature described by a 
temporal signature. The phrase temporal signature is generally 
construed as the characteristic of a feature over time. Ts give an 
extremely simplified example, one might say that a deciduous tree 
has the temporal characteristic of being light green in spring, 
richly deep green in summer, yellow in fall, and brown in winter. 

Temporal studies of film imagery were performed by creating 
color enhancements on film viewers using scenes of different dates 
on different channels. Temporal studies of digital tape data were 
performed by analyzing combined data sets of the same scenes at 
different dates (e.g., a composed data set of eight channels 
derived from registering two data sets of the same scene on two 
different dates). The capability existed on the Univac 1108 com- 
puters and on the ERIPS to register digital data sets; thus, com- 
posite temporal data sets were compiled. 

Atmospheric corrections and study of absolute reflectance . - 
The effects of the varying Sun angle and atmospheric conditions 
were expected to be significant and were known to change the 
detected radiance levels of various features. The removal or cor- 
rection of these effects was attempted to investigate the unique- 
ness of the spectral characteri'stics of the various features. 
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1 . 5 ACCOMPLI S HMENTS 

Significant improvements were made in existing data— processing 
and analysis techniques. Furthermore, the analysis techniques 
that were applied to the ERTS-1 data and to corollary information 
provided an appropriate framework for the systematic evaluation of 
important elements of a prototype multispectral data-processing 
and analysis system. 


Most of the techniques used in processing and interpreting 
ERTS-1 data were refined or developed in the context of resource 
application or use. Because most of these techniques must be 
closely tailored to operational requirements, participation by 
user organizations would have permitted them to become more knowl- 
edgeable in the development and use of techniques appropriate for 
their application. 


An outstanding example of this operational philosophy was the 
development of a computer-aided procedure to successfully detect 
and locate surface water using ERTS-1 data. The work was initi- 
ated at JSC in support of the State of Texas in its cooperative 
effort with the U.S. Army Corps of Engineers for an inventory of 
water impoundments pursuant to Federal legislation. This develop- 
ment of an accurate and efficient method for locating surface 
water using the ERTS— 1 data evolved from the investigations per- 
formed by the analysis teams. 

The procedure was evaluated for a test site area in east- 
central Texas and resulted in 100-percent identification of all 
lakes with a surface area of at least 4 square hectometers and 
location of those lakes within an accuracy range of 300 meters 
(1000 feet). This computer-aided procedure produces line printer 
output at any desired scale with symbols depicting the surface 
water. The printout is precisely registered to the base maps. 
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The procedures, with the data products, were provided to the Texas 
Water Rights Commission and Texas Water Development Board and will 
be evaluated with other techniques in Howard, Montgomery, and 
Washington Counties in Texas. This example is a specific instance 
of the transfer of remote-sensing technology to the user community. 

Using conventional interpretation and computer-aided process- 
ing, the feasibility of extracting useful resource information 
from satellite data was demonstrated as a result of this project. 
However, further refinements in the use of temporal data and 
radiometrically calibrated data, together with improved precision 
in compensating for atmospheric effects within the data, will 
increase the degree of accuracy to which selected surface features 
of applications interest can be detected and identified. 

The following is a summary of some of the significant find- 
ings that were developed in the course of this investigation. 

They are representative of the type of accomplishments that led to 
the formulation of the conclusions described in the next section 
of this document. 

1 . The usefulness of temporal data for detecting range fea- 
tures was demonstrated. Temporal analyses enhanced the discrim- 
ination between different classes; this conclusion conforms with 
mathematical theories . 

2. Level I land use signatures were extended over an area of 
4660 square kilometers {1800 square statute miles) with encourag- 
ing results. Data over the entire scene were sampled to develop 
classification training statistics instead of establishing train- 
ing fields throughout the scene. This approach holds promise for 
classifying Level I land use over very large areas for which 
ground truth is not available. 

3. Nonsupervised classification yielded information that 
could be related to various types of residential areas in part of 
Houston. 
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4. Computer-aided classification of a portion of the Sam 
Houston National Forest permitted identification of 14 classes of 
forest land conditions. 

5. Level II classification was satisfactory for separating 
wooded rangeland, nonwooded rangeland, cropland, and water. A 
Level III classification was possible in separating marshhay cord- 
grass and gulf cordgrass. 

6. The separation of crop types (Level T T I; i.e., small 
grains, truck farm crops, grasses, summer fallow) was accomplished 
by conventional interpretation as well as by computer-aided clas- 
sification to a 95-percent classification accuracy. Further 
breakdown to crop species (Level IV) was accomplished by computer 
methods only, but the accuracy was reduced to 70 to 80 percent 
fcr single-date observations. 

7. High spectral contrast between adjacent features permit- 
ted their distinction between each other. On the other hand, a 
"blooming effect" occurred in which the brighter element seemed to 
occupy a greater spatial extent than it should. This effect con- 
tributed to errors in areal measurements. 

8. The smallest field detected using conventional image 
interpretation was 6.5 square hectometers (16 acres). 

9. A technique was developed to eliminate the striping 
(banding) observed on all four channels of the ERTS-1 data and 
attributed to the residual differences (after GSFC calibration) in 
detector output. The striping is particularly evident when the 
scene is a large, homogeneous area such as a lake and is suffi- 
cient to affect sensitive, spectral-clustering algorithms. Use of 
the technique removes detector-to-detector differences while 
leaving the overall average value of the data unaltered. 

10. An extensive phytoplankton bloom was detected in the 
ERTS-1 band 6 imagery obtained over Galveston Bay, Texas, 

February 24, 1973. 



11. Changes in water coverage due to rice farming practices, 
flooding, and tidal-level changes were detected by temporal analy- 
sis of ERTS-1 data collected over the Trinity River delta and 
Trinity Bay area in Texas. 

12. Coastal features exhibited unique and repeatable reflec- 
tance signatures as obtained from ERTS-1 MSS data. These signa- 
tures or characteristic reflectance properties of coastal features 
were obtained by the application of a computer routine, which 

was developed during the investigation. 

13. Temporal signature extension was accomplished on large, 
relatively homogeneous features such as water bodies. Sun-angle 
changes insignificantly degraded temporal signatures with data 
sets separated by more than 36 days. Changes in turbidity caused 
the greatest difficulty in extending signatures of freshwater 
over time and distance. 


1.6 CONCLUSIONS 

Detailed conclusions are contained in the seven individual 
reports prepared by the investigation teams. The major conclu- 
sions drawn from these seven investigative areas are contained in 
this compendium in summary form. The highlights are as follows. 

1. Data of the type provided by the ERTS-1 multispectral 
scanner will be of significant value for conducting extensive 
regional inventories and surveys in the areas of land use, crop- 
land and range/fores t resources, and for determining certain 
coastal and estuarine resources and conditions such as water 
turbidity, marshland, forests, ancl phytoplankton blooms. 

2. Classification performance for various features using 
ERTS-1 data both by conventional image interpretation and by 
computer-aided methods compared favorably with classification 
performance historically achieved with high-altitude-aircraft 
data. 
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3. These investigations confirmed the assumptions that large 
crop and land resources inventories, as stated in USGS Circular 671 
for Level I land use survey, can be conducted with ERTS-1 data 
using relatively little ground truth. 

4. Field width, relative contrast, and orientation were 
important characteristics in accurately detecting the boundaries 
of individual fields. 

a. Narrow fields oriented east-west (parallel to the scan 
lines) were more difficult to detect and measure than fields 
oriented north-south. 

b. Fields less than 50 meters wide were not consistently 
detectable. 

5. Low computer-aided classification accuracies occurred 
when a spectrally complex urban scene was classified with exten- 
sive nonurban areas containing spectrally homogeneous features. 

To increase classification accuracies of certain urban land use 
categories, it was necessary to develop separate computer inputs 
and regroup some spectrally similar clusters. Even so, classifi- 
cation accuracies of urban landscapes did not approach the accura- 
cies achieved in classifying the land use categories containing 
more homogeneous features (agriculture, forest, water, etc.) 

6. Some categories of land use appeared to be best deter- 
mined by classifying at a detailed level and aggregating to more 
general levels. 

7. Sun-angle effects were significant in distinguishing an 
absolute reflectance and in class discrimination. Temporal analy- 
sis of ERTS-1 data sets demonstrated a high degree of signature 
overlap in the various spectral classes. Removal of Sun-angle 
effects resulted in good discrimination of the same classes. 
Atmospheric effect may not have a pronounced impact on spectral 
discrimination of classes in a relative sense, but must be removed 
to determine absolute reflectance. 
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8. Temporal analysis was shown to provide better discrimina- 
tion between certain features. 

9. The computer-aided classification techniques were gener- 
ally more accurate than conventional image-interpretation 
techniques for homogeneous features that could be differentiated 
by spectral contrasts. In separating urban and nonurfcan features, 
conventional image-interpretation methods proved advantageous for 
delineating small, heterogeneous features (urban, linear patterns, 
etc. ) . The choice of techniques may depend on the degree of 
classification accuracy acceptable, the computer time available, 
the skills of the analyst, and the quantity to be processed. 


2 • 0 THE ERTS-1 AGRICULTURAL ANALYSIS 
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2 . 1 OBJECTIVES 

The general objective of this investigation was to evaluate 
how well features of agricultural importance could be detected, 
identified, and located and their areal extent measured using 
ERTS-1 data. This general objective included the following 
specific objectives, which are listed in order of priority. 

1. To separate agricultural areas from nonagricultural areas 

2. To separate cropland from noncropland within the agri- 
cultural areas 

3. To determine whether the existence of different crop 

types (e.g., rowcrops and small grains) could be detected within 
the cropland 

4. To determine whether the existence of different crop 

species (e.g., wheat, barley, and oats) could be detected within 
the cropland 

5. For each major crop, to determine the size of the small- 
est field that could be detected’, identified, located, and 
measured 

6. To determine the effect of varying field shape and field 
size on the accuracy of crop classification and field measurement 

7. To determine the effect of the relative contrast of 
adjacent fields on boundary detection and crop classification 

A general classification hierarchy of the agricultural fea- 
tures is shown in figure 2-1. Additional objectives were to eval- 
uate the effect of atmospheric correction techniques on an 
agricultural analysis of ERTS-1 data and to become familiar with 
and to evaluate the capabilities of the JSC hardware, software, 
and procedures for processing and analyzing such data. 
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Figure 2-1.- Agricultural classification hierarchy. 
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2.2 ANALYTICAL APPROACH 
2.2.1 Data Processing 

The best ERTS-1 data set available from the 1972 crop year 
for each of six study areas was analyzed. The analysis included 
computer-aided and conventional image- interpretation techniques. 
During the analysis, one study area was selected for a temporal 
analysis. 

Area measurements of representative fields were made in the 
study areas using a mechanical-electronic planimeter to manually 
trace the field borders. Additive color-enhanced ERTS-1 imagery 
of approximately 1:130,000 scale was used for the measurements. 
Aircraft photography and other ground-truth data were used as a 
basis for comparison and for precise instrument calibration. 

The Universal Transverse Mercator (UTM) coordinate system was 
selected for tract location. The UTM grids with 1000-meter incre- 
ments were prepared for enhanced ERTS-1 imagery of approximately 
1:130,000 scale. The grids were constructed with a programed 
XY-plotter, which was set to scale by comparing the measured dis- 
tances between specific featv.res on the imagery to USGS topo- 
graphic maps. 


2.2.2 Study Areas 

Locations of the following six study areas within the United 
States are shown in figure 2-2. 

Hill County, Montana .- Hill County, Montana, was selected for 
study because of the unique field patterns that resulted from the 
practice of "strip-fallow" farming. The result is a series of 
long, narrow, rectangular fields 91 to 107 meters wide by 798 to 
1586 meters long. Wheat is the major crop in the area. However, 
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Figure 2-2.- Six agricultural study areas within the United States 
































barley, other small grains# alfalfa hay# and some cultivated 
pastures are also found. Data from August 7, 1972, were analyzed. 

Imperial County, California .- Imperial County# California, was 
selected because the yearlong growing season ensured a variety of 
crops in various stages of maturity (e.g., alfalfa hay, barley, 
sugar beets, cotton, grain sorghum, and winter vegetables). Most 
of the fields were well defined because they are bounded by roads, 
irrigation canals, or drainage canals. Data from November 6, 1972 
were selected for detailed analysis. 


Hardin County, Iowa .- The major crops of Hardin County, Iowa, 
a typical corn-belt farming area, are corn, soybeans, and oats. 

The earliest usable data pass occurred on August 13, 1972, with 
limited cloud cover in the north and south. The oats had been 
harvested by July, and a strong cover of secondary growth appeared 
mixed with some grasses and legumes. In most cases, the large 
blocks of corn represented a combination of several smaller fields 
with boundaries not visible from the ERTS-1 altitude. 

Holt County, Nebraska .- The Holt County, Nebraska, study site 
was unique because of its center-pivot irrigation systems. These 
contrasting targets of 53 square hectometers (130 acres) and a 
limited variety of crops made Holt County one of the easier sites 
to analyze. The test site contained approximately 30 fields of 
field corn, the major crop. In addition, there were three pop- 
corn fields, one sunflower field, several fields of alfalfa and 
fields of grass at various stages of growth, two fallow fields, 
and one sorghum field. Data sets from July 30 and August 16, 1972 
were used for analysis. 

Butte County, California .- Butte County, California, was 
selected primarily because rice was a major crop. Other crops in 
this county were fruit, nuts, and small grains. The fields were 


2-6 

large and irregularly shaped with medium contrast. Data from 
September 19, 1972, were selected for analysis. 


Worth County, Georgia .- The Worth County, Georgia, test site 
had approximately 30 percent of the 31-square-kilometer (12 square 
statute mile) study area in woods. The cultivated fields in the 
site were generally small, ranging in size from less than 
0.4 square hectometer (1 acre) to 40 square hectometers 
(100 acres). The majority of the fields were smaller than 
12 square hectometers (30 acres). Field boundaries followed 
natural contours, woods, and drainage patterns; thus, most of the 
fields were nonrectangular . Worth County was included in the 
analysis because it contained small and irregularly shaped fields 
of cotton, tobacco, and peanuts. The September 8, 1972, ERTS-1 
MSS data were selected for analysis. This set of data was the 
earliest available in which the study area was sufficiently clear 
of clouds. These data were not optimum because it was late in the 
season and many of the crops were already at the mature or harvest 
stage. 


2.3 RESULTS 

The results of the investigation showed that ERTS-1 data 
generally could be used for crop classification, location, and 
area measurement. Although applications procedures must be devel- 
oped, good classifications were achieved from the digital data 
processed by using clustering programs and the LARSYS program. 

Spatial information contained in the imagery was primarily 
used at classification Level I (agricultural compared to nonagri- 
cultural) . Generally, nonagricultural areas in the sites analyzed 
were irregular in shape, showed as line features, or had a reflec- 
tance characteristic similar to water or soil. Conversely, agri- 
cultural areas tended to exhibit regular shapes and high infrared 
(IR) reflectance if crops were growing. 


2-7 

Cropland was distinguishable from noncropland because virtu- 
ally all cropland either consisted of regularly shaped fields or, 
if irregular, the boundaries were well defined and obviously man- 
made. The presence of other clues, such as canals, location with 
respect to other fields, and similarity to known cropland, was 
also useful. The ground-truth maps of all the counties show the 
distinct field boundaries. 

The clustering routine 1S0CLS was useful for a rough Level I 
and Level II qualitative classification; that is, for the separa- 
tion of agriculture from nonagriculture and cropland from non- 
cropland. However, crop types or crop species were not 
satisfactorily separated in Worth County using this technique. 

The smaller fields blurred into a hodgepodge that was not identi- 
fiable. Field sizes of 10 to 12 square hectometers (25 to 
30 acres) were the minimum for computer crop identification. In 
this study area, approximately 30 of the total of 475 fields (less 
than 7 percent) were larger than 12 square hectometers (30 acres). 
If these fields had represented a cross section of the crops 
present, and if they had been at a stage of crop development 
making them suitable for use as training fields, fields as small 
as approximately 4 square hectometers (10 acres) might have been 
classified. Because of the inconclusive results of the classifi- 
cation, further analysis of the small fields of Worth County will 
not be discussed in this report. 

The most significant results for Levels III and IV, as shown 
in the classification hierarchy (fig. 2-1), are discussed in the 
following sections. These levels will be treated together under 
conventional image interpretation and computer classification. 

2.3.1 Conventional Image Interpretation 

Conventional interpretation of the enhanced imagery revealed 

. , that the relative contrast of adjacent features and their 

I t 

i 


geographical orientation were important in ascertaining their de- 
tectability. Linear features (such as long, narrow fields in Hill 
County) parallel to the scan lines on an MSS image were difficult 
to detect and define. The data for Hill County were acquired 
August 7, 1972, late in the growing season when most of the fields 
had already been harvested, and little contrast between the culti- 
vated fields and the fallow fields remained. Fields 90 meters 
(300 feet) wide that were oriented north-south were detectable. 

The narrowest field detected with an east— west orientation was 
135 meters (450 feet) wide. Table 2-1 includes the results of the 
crop identification and the percentage accuracy. 

TABLE 2-1.- HILL COUNTY CONVENTIONAL IMAGE INTERPRETATION 



Number of fields 


Summer fallow 


Barley 

Winter wheat 


Spring wheat 

Crested wheat 
grass 

Oats 

Totals 


Training Test Test fields Accuracy, 
1 1 identified 1 


percent 



Accuracy of field identification equals number of test 
fields identified divided by total number of test fields. 













2-9 


On the Imperial County test site, crops having similar spec- 
tral reflectance characteristics were arranged together into 
classes on the ACVP false-color IR enhancement (fig. 2-3). The 
smallest field easily identifiable was 6.5 square hectometers 
(16 acres) wide, rectangular, bright, and bordered by contrasting 
fields. A dividing border, such as a dirt road at least 40 meters 
(132 feet) wide, was required to detect the boundary between 
fields of the same contrast. Examination of figure 2-3 indicates 
a good breakdown at classification Level II (cropland compared to 
noncropland) , but further conventional interpretation was almost 
impossible. Using various training fields, correlation of a 
particular color to a single crop species was attempted. Again, 
the results quickly indicated that crop classification in this 
manner was almost totally ambiguous. The colors in figure 2-3 
were found to be an indicator of the density of vegetative cover, 
and these results appear in table 2-II. 

TABLE 2-II . - IMPERIAL COUNTY CONVENTIONAL IMAGE INTERPRETATION 


Color 

Crop species and height, cm 

Red 

Alfalfa 36, sugar beets 30, sudan grass 24 

Moderately red 

Alfalfa 21, 15, 9; sugar beets 30; 
asparagus 61 

Dull red 

Alfalfa 15, 9; sugar beets 9; carrots 21; 
sudan grass 21 

Pink 

Melon 15 

Magenta 

Sugar beets 27, 12 

White 

Bare soil; dry, plowed, or recently planted 

Black 

Bare soil; undisturbed with little or no 
ground cover 

Gray 

Sugar beets 9; grain sorghum stubble 61; 
bare soil; melons 15 
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The crops in Hardin County were identified in those fields 
that could be detected. The best color composite generated from a 
CCT at JSC (fig. 2-4) was used for crop identification. The ERTS 
band 6 was on the green gun, band 5 on the red gun, and band 4 on 
the blue gun. There is a cloud just south of the study area, and 
its shadow obscures the northern edge. Table 2-III shows the 
colors that were grouped for each crop type. The crop types were 
identified on the basis of color. These three crops (corn, soy- 
beans, and oats) and bare soil make up 85 percent of the study 
area. 

TABLE 2-III. - HARDIN COUNTY CONVENTIONAL IMAGE INTERPRETATION 




Number of fields 

Crop 

Color 

Total 

Training 

Test 

Test 

fields 

identified 

Accuracy, 

percent 3 

Corn 

Dark 

green 

89 

12 

77 

73 

95 

Soybeans 

Bright 

green 

75 

10 

65 

61 

94 

Oats 

Lighter 
red and 
pink 

62 

14 

48 

41 

85 

Bare soil 

Dark red 

3 

_! 

2 

2 

100 

Totals 


229 

37 

192 

177 



cl 

Accuracy of field identification equals number of test fields 
identified divided by total number of test fields. 


m Corn ■ Bare soil 
□ Soybeans lj Cloud 
mz\ Oats HZ) Cloud shadow 


Figure 2-4.- Hardin County conventional image interpretation of JSC color composite 

(NASA S-73-281 41 ) . 
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2.3.2 Computer Classification 

Because of the lack of ground data, atmospheric correction 
techniques were not used in the computer-aided analysis. Photom- 
eters, which would have measured the effect of the atmosphere on 
the solar energy reaching the Earth, were not available for the 
test sites during the ERTS-1 overpasses. 

Hill County .- All the major crops grown on the Hill County 
site had either been harvested or were ready for harvest. In 
addition to the advanced stage of maturity , the narrow width of a 
majority of the alternating strips made it very difficult to 
select good training and test fields of sufficient size to achieve 
accurate classification. 

An analysis of the results from the Hill County site in 
August indicated that fields smaller than 24 square hectometers 
(60 acres) that were longer than 800 meters (0.5 statute mile) 
could not be properly defined for use as training or test fields 
on the interactive classification system. These narrow fields had 
a higher percentage of picture elements (pixels) comprising the 
boundary and therefore might not represent the reflected data from 
that field. Many fields in the area were 1856 meters (1 nautical 
mile) long, which meant they must be larger than 40 square hecto- 
meters (100 acres) to be wide enough for use as a dependable 
training or test field. Very few of the fields could meet these 
criteria. In contrast, fields that were as small as 8 square 
hectometers (20 acres) made dependable training or test fields 
provided the strips were less than 400 meters long and adjacent to 
summer fallow strips. The high relative contrast between the 
small grain crops and summer fallow strips provided an excellent 
boundary discrimination until the crop was harvested, at which 
time the contrast was greatly reduced. The accuracy of classifi- 
cation in Hill County for the training and test fields is shown 
in table 2-IV. 


2-14 


TABLE 2-IV. - HILL COUNTY CLASSIFICATION PERFORMANCE SUMMARY 
(a) Barley, spring wheat, and winter wheat 


Class 

No. 

samples 

Barley 

Spring 

wheat 

Winter 

wheat 

Fallow 

Sod 

Correct 

classification, 

percent 


Training field 

Spring wheat 


33 

93 

9 

7 

1 

65.0 

Winter wheat 

184 

34 

8 

135 

4 

3 

73.4 

Barley 

147 

129 

9 

4 

4 

1 

87.8 

Fallow 

188 

7 

3 

3 

164 

11 

87.2 

Sod 

173 

2 

3 

3 J 

13 

155 

89.6 




Test 

field 




Spring wheat 

52 

11 

35 

6 

0 

0 

67.3 

Winter wheat 

99 

18 

49 

27 

5 

0 

27.3 


{b) Barley and wheat 


Class 

No. 

samples 

Barley 

Wheat 

Fallow 

Sod 

Correct 

classification, 

percent 




Training 

field 



Wheat 

327 

67 

245 

11 

4 

74.9 

Barley 

147 

129 

13 

4 

1 

87.8 

Fallow 

188 

7 

6 

164 

11 

87.2 

Sod 

173 

2 

6 

13 

155 

89.6 

Test field 

Wheat 

151 

29 

117 

5 

0 

77.5 


(c) Small grains 


Class 

No. 

samples 

Small grains 

Fallow 

Sod 

Correct 

classification, 

percent 



Training 

field 



Small grains 

MM 

454 

15 

5 

95.7 

Fallow 


13 

164 

11 

87.2 

Sod 

n 

8 

13 

155 

89.6 


Test field 


Small grains 

151 

146 

5 

0 

96.7 













































t 


t 
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The classification results are best shown by the comparison 
of the ground-truth map for the 31 -square-kilometer (12 square 
statute mile) test site (fig. 2-5) with the recognition map 
developed (fig. 2-6). Crop identification accuracy was computed 
using the following formula: percent accuracy of field identifi- 

cation equals number of samples (pixels) per feature (test or 
training field) identified divided by the total samples (pixels) 
per field (test or training) multiplied by 100. Oats were grown 
in the Hill County site; however, no classification attempt was 
made because the crop had been harvested at the time of the over- 
pass and no oat field was large enough for use as a training 
field. 

Imperial County .- Using the NSCLAS clustering routine and 
ground-truth information, eight different crop types and species 
were selected for classification training and testing in Imperial 
County. Although corresponding test fields were not available, 
additional training fields were selected so that the classifica- 
tion might be a fairly accurate representation of the actual 
ground truth as shown in figure 2-7. Table 2-V contains the per- 
centage accuracy figures for the training and test fields outlined 
in yellow on figure 2-8. 


The low percentage of classification accuracy for beets may 
be attributed to their being at a very early growth state. More 
ground was visible at this stage, and beets were misclassif ied 
as bare soil and sorghum. The low percentage of classification 
for bare soil was due to the similarity in spectral response of 
sorghum and bare soil. 



□ Sod 

LJ Summer fallow 
E Summer wheat 
E3 Winter wheat 
O Barley 

■ Oats 

□ Crested wheat grass 

■ Water, building sites 


Section identification 


13 

18 

17 

16 

15 

14 

24 

19 

20 

21 

22 

23 


Roads — 


Figure 2-5.- Hill County 1972 annual ground truth (NASA S-73-23401) 
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Fi<jure 2-6.- Hill County computer classification (NASA 
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Mr? • 


[_J Pasture 
f [ Beets 
H Melons 
LJ Hemp 

Q Grain sorghum 


■ Alfalfa 
H Bare soil 

0| Wheat stubble 

■ Sudan grass 
Hi Sugar beets 


| Seed onions (planted) 
H Asparagus 
Q Sesame 
[_| Citrus 
■ Cotton 


[J Lettuce 
£|§ Carrot 
Q Weeds 
Tillage 


Figure 2-7.- Imperial County 1972 annual ground truth (NASA S-73-25512) 






TABLE 2-V. - IMPERIAL COUNTY CLASSIFICATION PERFORMANCE SUMMARY 
(Nov. 6, 1972] 


designated |£rrSct Carrots bo!!®! ,011*2 Beets Alfalfa As P aragUS Cotton Sorghum Melons 

group r 


raining field performance 



Carrots 

32 

7.. 9 

Bare soil 1 

32 

65.6 

Bare soil 2 

70 

90.0 

Beets 1 

20 

95.0 

Beets 2 

25 

0.0 

Beets 1 

45 

42.2 

and 2 



Alfalfa 

24 

95.8 

AsparaguB 1 

80 

95.0 

Asparagus 2 

56 

55.4 

Asparagus 1 

136 

78.7 

and 2 



Cotton 1 

27 

63.0 

Cotton 2 

48 

100 

Cotton 1 

75 

86.7 

and 2 



Sudan grass 

24 

91.7 

Totals 

* )» 

(b) 


Teat field performance 



Overall performance: 264/286 » 92.3 percent 

Overall performance: 343/438 » 78.3 percent 


































































CU Alfalfa 
H Asparagus 
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I I Sudan grass 
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d Bare soil 1 
I I Bare soli 2 
H Beets 


Imperial County computer classification 
1972, data (NASA S-73-28070) . 


November 6 
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Hardin County .- Corn was the dominant crop in Hardin County 
and was readily recognized in a supervised classification. A 
total of 76 cornfields was represented on the ground-truth tab- 
ulation (fig. 2-9). Of these, 73 were classified correctly. The 
three fields unaccounted for were small, individual fields, beyond 
the expected resolution limit. The larger fields were fairly well 
defined with little extraneous data. The thin cloud to the south 
of the study area cast a weak shadow along the northern edge of 
the area and caused other crop types in this area to be grouped 
with the corn. One corn test field was partly shadowed. Removing 
the data for the shadowed field increased the corn test results 
to higher than 80 percent. Table 2-VI contains a breakdown of 
eight classes. 

The classification results were reasonably good an a pixel- 
by-pixel basis and indicated that corn, soybeans, and oats could 
be separated at this maturity stage (table 2-VI) . A review of the 
input showed that in almost all cases of low classification accu- 
racy , one of three things happened. The most frequent occurrence 
was improperly defined training field boundaries overlapping 
another crop type. Because this area had many small fields 
(fig. 2-9), this problem was expected,. Judicious screening of the 
training and test fields as shown in pink on figure 2-10 overcame 
problems with smaller fields. The second problem area was that 
all oats had been harvested nearly a month before the overpass and 
were in various stages of regrowth or in a bare-soil condition. 
Despite this condition, the classification was fairly successful, 
except for an inability to distinguish some oats from county and 
state roads. This difficulty was probably attributable to the 
fact that farm roads in Iowa are normally bordered by wide, shal- 
low ditches of mowed bluegrass, which appeared spectrally similar 
to harvested oats. The third problem was presented by the cloud 
shadow that was cast over the northwest corner of the study area. 
No special attempt was made to develop unique signatures for this 


area. 



Brome 


Feed lots 


Foxtail 


Section identification 


Soybean 


Sorghum-sudan 


Plowed 


Clover 


Alfalfa 


Farmsteads- roads 


Alfalfa-oats-clover 


Pasture 


Alfalfa-brome 


Brush -weeds 


Hardin County 1972 annual ground truth (NASA S-73-25510) 


17 

16 

15 

14 

13 

18 

20 

21 

22 

23 

24 
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TABLE 2-VI .- HARDIN COUNTY CLASSIFICATION PERFORMANCE SUMMARY 

[Aug. 13, 1972] 


Field 

No. 

samples 






Bare 

soil 




designated 

group 

Corn 

Soybeans 

Oats 

Sorghum 

Sudan 

Cloud 

Shadow 

Accuracy , 
percent 


Training field performance 


" 

Corn 

419 

380 

1 

8 

21 

0 

2 

0 

7 

90.7 

Soybeans 

267 

0 

220 

1 1 

17 

15 

0 

0 

4 

82.4 

Oats 

106 

2 

2 

86 

4 

1 

10 

0 

1 

81 .1 

Sorghum 

56 

0 

3 

1 

52 

0 

0 

0 

0 

92.9 

Sudan grass 

47 

0 

4 

0 

2 

41 

0 

0 

0 

87.2 

Bare soil 

58 

3 

0 

6 

0 

0 

49 

0 

0 

84.5 

Cloud 

491 

0 

0 

2 

0 

0 

0 

489 

0 

99.6 

Shadow 

1 31 

1 

3 

0 

1 

0 

0 

0 

126 

96.2 

Test field performance 

Corn 

309 

223 

6 

14 

14 

0 

9 

0 

a 43 

72.2 

Soybeans 

141 

1 

128 

0 

5 

5 

0 

0 

2 

90.8 

Oats 

85 

1 

5 

64 

4 

2 

8 

1 

0 

75.3 

Cloud 

365 

0 

0 

1 

0 

0 

0 

364 

0 

99.7 

Shadow 

98 

5 

0 

0 

3 

0 

5 

0 

85 

86.7 


Corn test field partly shadowed; removal of these data increased corn test results to more than 
80 percent. 










































E'iyure 2-10.- Hardin County supervised computer classification, Auqust 13, 1972 
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Holt County .- The fields in Holt County were analyzed in two 
phases (fig. 2-11). In the first phase, four channels of ERTS-1 
MSS data were analyzed for a single pass on July 30, 1972 
(fig. 2-12 (a)). In the second phase, the August 16, 1972, imagery 
was correlated and registered with the earlier imagery to form a 
single data set. Data from the second phase were analyzed using 
both a 0.5 and a 0.2 threshold. The 0.2 threshold was the best 
(fig. 2—12 (b) ) . The results of the first and second phases of 
analyses are shown in table 2-VII. These results are summarized 
by crop species as follows. 

Field corn or popcorn: The classification accuracy for the 

identification of known test cornfields was approximately 97 to 
98 percent for both data sets. The single-pass data set did not 
contain the information necessary to separate popcorn from field 
corn. An attempt to distinguish field corn from popcorn using a 
two-pass data set was highly successful. This differentiation was 
possible because the reflectance of popcorn fields in channels 3 
and 4 of the ERTS-1 MSS in the August 16 pass was higher than that 
of field corn. Such a spectral response may well be due to some 
peculiarity in the development of popcorn that is different from 
that of field corn. For example, popcorn matures faster than 
field corn. 

Sunflowers: Sunflowers were separable from corn because cf 

higher reflectance in the third channel of the single-pass data 
sets of July 30, 1972. At this stage, the sunflowers were bloom- 
ing. Sunf lowers were not separable from corn, particularly pop- 
corn, using the August 16 data only. 


ISJ 

I 

ro 

o 



1 Corn 

■ 

Alfalfa 

□ 

Sudan 

[ i Popcorn 

□ 

Pasture 

m 

Grain sorghum 

(3 Grass 

□ 
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— 

Section line 

H Fallow 

■ 

Christmas trees 



[_] Sunflowers 

H 

Weeds 




Section identification 


18 

17 

16 

15 

14 

13 

19 

20 

21 

22 

23 

24 


Figure 2-11.- Holt County annual ground truth 





Field 

designated 

group 



TABLE 2-VII . - 

HOLT COUNTY CLASSIFICATION PERFORMANCE 

: SUMMARY 



Data set 1 (single pass), 3 
threshold = 0.5 

Data set 2 (two pass) . 
threshold = 0.5 

1 ] 1 1 

Data set 2 (two pass) 
threshold = 0.2 

-i r — : 

b 

t 

“I 


Total 

samples 


Thresh- 

Misclassif ied old 


Percent Total 
correct samples 


Thresh- 

Misclassif ied old 


Percent Total 
correct samples 


Misclassif ied old 
out 


Percent 

correct 


Field corn 
or 

popcorn 
Field corn 
Popcorn 


Sunflowers 

51 

Alfalfa 

34 

Alfalfa 

and 

grass 

25 

Grass 

96 

Pasture 

232 

Fallow 
(bare 
soi 1 ) 

34 

Grain 

sorghum 

— - 


Field 

corn or 
popcorn 

Field 

corn 

Popcorn 

Grass 


4 

99.0 

1038 




768 

— 

-- 

270 

0 

98.0 

49 

0 

85.3 

29 

0 

6 B . 0 

29 

1 

95.8 

97 

1 

91.8 

21 1 

0 

100.0 

44 

— 

— 

— 


Training field performance 


JO ! 7 Ilf 


Test field performance 


0 

97.5 

47B 

— 

— 

428 


— 

50 

0 

100.0 



50 


98.8 

979 

94.7 

745 

C 95 . 3 

234 

98.0 

50 

91 .2 

34 

96.6 

18 

94.8 

1 09 

98.2 

219 

100.0 

50 


14 


a Four -channel classification, July 30, 1972. 
b Seven-cher,r.el classification, August 16 and July 30, 1972. 

c Value does not reflect that one of seven known popcorn fields was misclassif ied as field corn. 
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Corn 

Hay grass 

Pasture 

Sunflower 


I — J Alfalfa 

I — I Alfalfa and grass 

■§ Bare soil 


(a) Supervised computer classification map of six classes using 

July 30, 1972, data. 


I 


Figure 2-12.- Holt County computer classification (NASA S-73-1506). 
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(b) Temporal composite using July 30 and August 16, 1972, data. 


Figure 2-12.- Concluded (NASA S-73-1506). 



2-30 


Alfalfa: The study site contained only a few small fields of 

alfalfa, some of which had recently been harvested at the time of 
both ERTS-1 passes. Alfalfa fields were classified as either pas- 
ture or grass in the second pass. Specifically, for the data of 
the first ERTS-1 pass, only two fields of usable size were identi- 
fied at or near the blooming stage during the second pass. These 
two fields were used as training fields for both data sets. To- 
gether, they comprised a total of approximately 30 data points. 

The accuracy of classification for these fields was 85 percent 
for the single-pass data set and 90 percent for the two-pass data 
set. 


Pasture and grass: Pasture and grass were easily separable 

from the other classes of the study site and were easily differ- 
entiated. The pasture was divided into three subclasses, one of 
which was bromegrass. The distribution of the bromegrass data 
overlapped the other types of grasses in the study site, in some 
cases by as much as 50 percent. Bromegrass was therefore not 
separable from other grasses in the area. 

Fallow land (bare soil): Bare soil has a unique spectral 

response and was easily separable from all types of vegetation. 
The two fallow fields in the study site were assigned to the same 
class, although the combined distribution of the data was defi- 
nitely biraodal. No difference between the fields was discernible 
in the classification maps, which were generated by assigning the 
two fallow fields to either the same or different subclasses. 

This similarity was true for the two-pass data set, even though 
some of the fallow fields seemed to be overgrown with vegetation 
during the second pass. Positive recognition of fallow fields 
may require that one of the two passes be made when the field is 
primarily bare soil. 
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Butte County .- Butte County displays a large variety of crop 
types. Analysis of the data from the ERTS-1 pass of September 19, 
1972, revealed that corn, rice, and fallow land were the most 
separable of the 17 classes used for training. Only 8 of these 
17 classes were selected for testing classification accuracy 
because o' the limited representation of the remaining classes in 
the test site (table 2-VIII) . 

Three classes had classification accuracies of less than 
70 percent: corn, native trees, and weeds. The low classifica- 

tion of corn could be partly attributed to the weeds and drowned 
spots in fields, which caused corn to be confused with pasture. 

Ten percent of the weedy corn (corn) was classified as good corn 
(corn 1). When the misclassif ications into pasture and corn 1 are 
disregarded, corn had a classification accuracy of 84 percent. 

The low classification of weeds resulted from confusion with 
alfalfa and trees (plum, almond) . The similarity of weeds and 
alfalfa was somewhat understandable. However, the confusion of 
weeds with trees can only lead to the conclusion that information 
was needed in the ground-truth reports to describe the density of 
the native trees. The five rice test fields had several distinct 
stages of growth, which accounted for the wide range (0 to 92 per- 
cent) of classification accuracy. The zero classification was in 
a test field where harvesting had begun, and the field was classi- 
fied as wheat stubble. The other test crops (plums, corn, beans, 
and almonds) had very low classification accuracies that were not 
readily explainable. 


TABLE 2-VtU.- BUTTE COUNTV CLASSIFICATION PERFORMANCE SUMMARY 



Overall per fornancus 834/1078 * 77.4 percent 
Overall performance! J9T/6U7 » 60.4 percent. 
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2.3.3 Area Measurement 

All area measurements were performed on imagery generated on 
the ACVP. The results were good with some qualifications. The 
Imperial County area measurement results are summarized in 
table 2-IX. Boundary definition was the key to accurate field 
measurement. Many fields were too small or too narrow to be 
clearly defined. Rectangular fields, especially long, narrow 
ones, that were oriented nearly parallel to the image scan lines 
were difficult to define. Accurate boundary definition was also 
difficult at high-contrast borders, such as those between a field 
with low spectral response and a field with high spectral response. 
This situation caused a blooming effect, in which the field with 
high response appeared larger and the field with low response ap- 
peared smaller. In Hill County, dimensions greater than 220 meters 
under optimum conditions of contrast and orientation were neces- 
sary for clear definition. 


2.3.4 Field Location 

Specific tracts were located by correlation to a UTM grid 
system. A technique was developed in which a UTM grid was over- 
laid on an ACVP enlarged image and matched to control points with 
known UTM coordinates on 1 : 62 , 500-scale maps. Slight variations 
in the scale of the enlarged image were compensated for by a 
family of grids, each varying slightly in scale. Matching the 
appropriate grid to the control points made possible the location 
of tracts in the 31 -square-kilometer (12 square statute mile) 
study areas to within approximately 200 meters. 

2.4 CONCLUSIONS 

The following conclusions can be stated. 

1 . A hierarchal approach was generally followed in this 
investigation. Agricultural areas were easily separated, and 
cropland within these areas was easily identified. These general 
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TABLE 2-1X.- IMPERIAL COUNTV FIELD AREA MEASUREMENTS 
[EHTS MSS imagery, Nov. 6, 1 9? 2 1 



r~ 


— 1 



M 

easured 

area 

Actual 

area 

P 

ercent 


Field 

designation 

Description 

arm/field - 
number 

. 2 
hm 

acres 

. 2 
hm 

acres 

error 


A 

mall, square field with 

/A 142-6 

13.2 

32.7 

13.7 

33.8 

-3.3 * 

ilfalfa; good growth and 1 
good border on two | 


strong vegetation re- 
sponse? two good borders 







sides 


and two low-contrast 
borders 

19/A! 60- 1 

13.1 

32.3 

15.1 

37.4 

-13.6 

lelons; poor border on J 

B 

Small, square field with 

three sides 

moderate vegetation; 







j 


strong vegetation re- 
sponse on two borders. 









low contrast on two 
borders 

49 1/A270- 1 


39.3 


36.9 

6.5 

Alfalfa; good growth and 

C 

Small, square field with 





borders j 


strong vegetation re- 
sponse; good borders on 









three sides 







Bare soil; poor west 
border 

D 

Very small, rectangular 
white field; two bor- 
ders poorly defined 

396/A436-4 

6.1 

15.0 

6. 3 

15.6 

-3.8 

Alfalfa: good growth and 


349/C511-4 


50.6 

19.3 

47.6 

6.3 

E 

Triangular field with 
strong vegetation; 
good borders 

9.8 

-32.1 

borders 

Melons; poor north border 





5B7/C6- 1 3 

24.1 

14,4 

35.5 

F 

Small, square field with 
moderate vegetation; 
good border on two 
sides, poor border two 
sides 

Alfalfa (7): good growth 







587/AUU-tl 

29.5 

73,0 

28.3 

70.0 

4 . 3 

G 

Rectangular, moderate size 
field with strong vege- 
tation response and 
good borders 

-1 , 2 

Dare soil; gpod borders 






587/A64-1 


72.1 

29.5 

73.0 

H a 

Rectangular, moderate size 
field of bare soil; good 
borders 

20.1 

3.3 

Alfalfa; weak growth and 





1 404/C51 5-4 

69.5 

27.2 

67.3 

I 

Rectangular, moderate si2e 


good borders ! 


field of low vegetation 








J a 

response; good borders 

Large, square fields cf 
bare soil; various bor- 
der contrasts 

1553/CS45-1 

59.3 

146.5 

60.1 

148.6 

-1 .4 

Bare soil; good borders I 
Light vegetation; diffi- 1 


745/C304-1 


36.4 

14.2 

35.2 

3.4 

K 

Small, square fields of 



cult borders i 


low vegetation re- 
sponse; various border 









contrasts 







Asparagus; good border ! 

1 L 

Large, square field of 

1136/C308- 

57.2 

141.3 

58.4 

144.4 

-2 1 


moderate vegetation 
response? good border 








N ♦ M 

Large, square fields of 
bare soil; various 
border contrasts 

877/C342- 1 
877/Cl 52- 1 

183.0 

452.2 

177.5 

438.6 

3. 1 

Bare soil; good border j 

Carrots; good to moderate! 


e 102/C65-1 

31.4 

77.5 

31.3 

77.4 

. 1 

S 

Rectangular, moderate siz 


borders 


field of low vegetation 
response; various 
border contrasts 







Asparagus, total; good 
borders 

P + R 

Large, square fields of 

1136/C30B- 

3 180.8 

446.7 

178.2 

440. 

1 . 4 


' moderate vegetation 

A26 3- 1 








response; moderate to 
poor borders 

A2B2-2 





— 

J > 


a Used for image scale calibration, scale = 1:130,000 
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separations were made primarily from spatial, rather than 
spectral, information gleaned from the ERTS-1 data. The separa- 
tion of crop types (i.e., small grains, truck farm crops, grasses, 
summer fallow) was accomplished by conventional interpretation as 
! well as by computer classification. 

2. Conventional image interpretation was valid for Level I 
and Level II classification with at least 95-percent classifica- 
tion accuracy. Further breakdown into Level III and Level IV was 
accomplished, but the accuracy was reduced considerably when an 
identification of individual species was attempted. 

3. For image-interpretation classification techniques, 
recombined false-color IR renditions were superior to all other 
additive false-color combinations that were attempted. Image 
enhancement and image interpretation revealed that colors related 
more to the density of vegetative cover than to crop type. 

I »• Boundary detection was found to be the key to field 

f detection and field measurement. Field width, relative contrast, 

t; 

I and orientation played an important part in the capability to 

j; accurately detect the boundaries of individual fields. In Hill 

(' County (where there are many long, narrow fields) , fields less 

- 

than 220 meters wide were not consistently detectable. Likewise, 
l fields that had an east-west orientation (parallel to the scan 

I lines) were much more difficult to detect and measure than those 

with a north-south orientation. Adjacent fields showing high 
relative contrast had boundaries that were very difficult to 

i 

define. This condition was caused by what is known as the bloom- 
j ing phenomenon. Likewise, adjacent fields showing low reflectance 

I were also difficult to define. The smallest easily detected 

! field in the Imperial County study area was 6.5 square hectometers 

‘ 1 

; (16 acres). On an enhanced image, the small field appeared white 

. i 

1 and rectangular, bordered by contrasting fields. 

at! 

i 

! 

i 

i 

f 

i 
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5. Computer classification of the ERTS-1 data was demonstra- 
ted to be very successful for many important crops. Clustering 
maps and ERIPS displays with baseline photography and maps of the 
subject were the best tools for selecting training and test fields, 
crop identification was accomplished to Levels III and IV in five 
of the test sites. Fields smaller than 12 square hectometers 
(30 acres) were not easily identifiable. Worth County, Georgia, 
had many such small fields, and the classification results were 
very poor. The long, narrow test fields of Hill County produced 
rather poor classification accuracies (less than 70 percent) . The 
smaller fields of Imperial County had an overall test-performance 
average of 78.7 percent. Hardin County, which has larger fields 
with rather poorly defined field boundaries, had a test field 
accuracy of 79.1 percent. Butte County, California, with fields 
having a variety of shapes and sizes, had an overall test field 
accuracy of 60.4 percent. Holt County (with large, well-defined 
fields) had a test field accuracy of 98 percent. 

6. Features on the surface could be located on ERTS-1 
imagery by correlation with a UTM grid system. System-corrected 
(bulk) data in conjunction with base maps were suitable for this 
procedure. For each study area, the features were located to ^ 
within approximately 200 meters by overlaying a grid on an AC VP 
enhanced image. This method is suitable for any test site 
covered by ERTS-1 for which large-scale base maps are available. 
Area measurements can be made of fields as small as 6.5 square 
hectometers (16 acres) under optimum conditions. In Imperial 
County, measurement error was in the ±3- to 4-percent range for 
fields having moderate to good border definition. In Hill County, 
the error was higher because of the long, narrow fields. 

7. A general conclusion can be reached that only the large 
(12 square hectometers or more) well-defined fields should be con- 
sidered for Level III and IV classification with the existing 
programs. However, with increased knowledge of the computer pro- 
grams and of parameters affecting the signatures and with improved 
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spectral information (regardless of source) , it is possible that 
Level III and IV classification can be obtained with reliable, 
repeatable accuracy. 
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3.0 THE ERTS-1 COASTAL/ESTUARINE ANALYSIS 

3.1 OBJECTIVES 

The specific objectives of this investigation were concerned 
with determining the performance of the ERTS-1 system in obtain- 
ing data capable of being used to map certain coastal features. 

To produce maps of a feature, one must be able to detect, iden- 
tify, and locate the feature and to determine its areal extent. 

The subject coastal features are Ixsted in figure 3-1 in a hier- 
archy designated Levels I, II, and III with each successive level 
corresponding to subclasses within a more general class. The 
hierarchy shown is not intended to be exhaustive but to show the 
major classes of features selected for this investigation. 

3.1.1 Classification Accuracy 

The first specific objective v/as to determine the accuracy 
with which specific features could be classified using ERTS-1 
data. Refer to tables 3-1 and 3-II. 

3.1.2 Feature Detection and Location 

The second specific objective was to determine the accuracy 
to which processed ERTS-1 data could be used to locate the posi- 
tion of detectable study features. 

3.1.3 Spectral Uniqueness and Temporal Ana: -.is 

The third specific objective was to determine the spectral 
uniqueness of detectable features in the ERTS-1 spectral bands as 
a function of time. This objective was concerned with the iden- 
tity of coastal features from the spectral and temporal properties 
of the feature as seen by ERTS-1 systems. 



Level I 


Level H 


Level m 



Figure 3-1.- Feature hierarchy for coastal investigation. 
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TABLE 3-1.- ACCURACY OF AREAL MEASUREMENT OF SALT MARSHES 


Marsh 

2 

Area, hm 

Accuracy 

, percent 

Aircraft 

data 

ISOCLS 

NSCLAS 

ISOCLS 

NSCLAS 

1 

351 

350 

353 

99.8 

99.4 

2 

77.4 

79.6 

76.6 

97.0 

99.0 

3 

60.5 

60.9 

67.1 

99.2 

89.0 


TABLE 3-II . - ACCURACY OF CLASSIFICATION OF 
BARRIER FLAT VEGETATED AREAS FOR 
GALVESTON ISLAND 
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3.1.4 Areal Measurement 

The fourth specific objective was to determine the accuracy 
to which processed ERTS-1 data could be used to measure the area 
of detected and identified features. 


3 . 2 ANALYTICAL APPROACH 

The investigators analyzed ERTS-1 data acquired in the area 
of Galveston Bay, Texas, including the surrounding wetlands, bar- 
rier islands, and coastal mainland areas. The basis for the 
accuracy determinations of coastal feature mapping was high- 

altitude photography from the NASA WB-57F and NC-130B aircraft as 
supported by field surveys. 

The selected coastal study area was viewed by ERTS-1 for 2 
consecutive days at approximately 10:22 local standard time every 
18 days, beginning August 28, 1972. The ERTS-1 data used for 
computer-aided and conventional analysis consisted of that 
acquired August 28 and 29, October 3, and November 26, 1972, over 
the Texas coastal study site. In addition, ERTS-1 data acquired 
January 19 and February 24, 1973, were used for conventional 
analysis. Since clouds block visible and near-infrared radiation, 
the team analyzed only the data acquired on predominantly cloud- 
free days. Of the ERTS-1 passes over the Texas coastal study 
site, 29 percent were on relatively cloudfree days. 


In addition to conventional image interpretation of aircraft 
and ERTS-1 images, data received from GSFC on computer-compatible 
tapes were reformatted and subjected to several automatic pattern 
recognition techniques and postprocessing procedures. The auto- 
matic pattern recognition techniques included both supervised 
(ERIP and LARSYS) and nonsupervised {ISOCLS and NSCLAS) classi- 
fication. After analysis and classification, the display options 
included cathode-ray tube, computer printout, and film trans- 
parencies. In subsequent postprocessing, an attempt was made to 


convert ERTS grayscale radiance to reflectance by using the 
ROTAR atmospheric correction algorithm. 

3.3 ANALYSIS RESULTS 

3.3.1 Water Features 

Water areas were discriminated easily in the ERTS-1 data 
through the use of a nonsupervised classification algorithm that 
defined natural spectral information units in the data using all 
four channels. From an analysis of ERTS-1 data, surface turbid- 
ity was assessed to be the property of water most strongly indi- 
cated by water-reflectance differences. The maximum solid matter 
(particulate or sediment) content in the water of the coastal 
study area ranged between 20 and 120 p/m. Ground-truth radiance 
measurements using an Exotech radiometer indicated a strong 
linear relationship between turbidity and radiance levels in the 
500- to 600- and 600- to 700-nanometer spectral bands, whereas 
- the 700- to 800- and 800- to 1100-nanometer bands suggest a 
curvilinear relationship. 

An example of the turbidity patterns in the ERTS-1 data is 
given in figure 3-2 in which land is red, wetland is yellow, and 
water turbidity zones are displayed by a sequence of colors from 
dark blue to medium blue, light blue, light green, and dark 
green to light brown. The least turbid water is dark blue, and 
the most turbid water is light brown. This display was processed 
on an EOD DAS using the results of an unsupervised classification 
analysis of the digital ERTS-1 data acquired November 26, 1972, 
over Galveston Island, Texas. 

One set of imagery acquired February 24, 1973, shows a 
phytoplankton bloom (reddish-brown algae, Exvella balteca) in 
Galveston Bay. Figure 3-3 is MSS band 4, and figure 3-4 is MSS 
band 6. Note the bright area in the western half of Galveston 


I 



land 

Wetland 

Least turbid water 

1 

Most turbid water 


EZ3 
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Figure 3-2.- Computer classification map of area surrounding 
Galveston Island, Texas, showing turbidity zones (ERTS-1 data 
acquired November 26, 1972) (NASA S-73-28182). 
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Bay in figure 3-4. This area is correspondingly dark in fig- 
ure 3-3. The Texas Parks and Wildlife Service determined that the 
bright area was a nontoxic bloom of the reddish-brown algae that 
began spreading from Seabrook, Texas, in late January to the 
imaged occurrence on February 24, 1973. Also, a suspected area 
of oil pollution in Galveston Bay was detected by the nonsuper- 
vised classification algorithm in the ERTS-1 data obtained 
August 29, 1972. This detection was manifested by the misclassi- 
fication of a known water area as salt marsh. 

3.3.2 Wetland and Dryland Features 

A number of wetland and dryland features of applications 
interest have been detected, identified, and mapped in the ERTS-1 
data through the use of a nonsupervised classification algorithm. 
Figure 3-5 is a display of the classification map of the Galveston 
data. A comparison of figures 3-5 and 3-6 is invited to determine 
the excellent agreement between the processed ERTS-1 data and the 
aircraft imagery acquired November 7, 1972 (fig. 3-6). 

3.3.3 Spectral Uniqueness 

The analysis of the ERTS-1 data also demonstrated that 
coastal features exhibit unique and repeatable reflectance signa- 
tures. From a plot of the ERTS-1 reflectances of all coastal 
features (all 3 months and all three sites) in MSS bands 5 and 7 
(fig. 3-7) , it can be seen that the signatures are unique, even 
with time. The partitions in figure 3-7 are valid only for the 
data in the figure and may change when more data are added. The 
closeness of signatures in some cases indicates problem areas for 
spectral classification. The corresponding plot of grayscale 
signatures (fig. 3-8) shows that grayscale signatures are not 
unique except for water and most sand areas. The primary correc- 
tion on the data was determined to be that of Sun angle. 
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3— 5.— Classification map of areas surrounding Galveston 
Island, Texas, showing wetland and dryland features. Nonsuper- 
vised classification of ERTS-1 data acquired November 6, 1972 
(NASA S-7 3-28 1 54 ) . 




gure 3-6.- Aircraft photograph of area surrou 
Island, Texas. Imagery acquired November 7, 
MB-57F aircraft at an altitude of 18 288 mete 
RC-8 metric camera, 2443 estar base Aerochrom 
(NASA S-73-28 1 0 1 ) . 
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3 ~ 7 *“ The ERTS ~} reflectance signatures of coastal features in MSS bands 5 
before corrections were made. (Data taken from all cluster analyses.) 
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Figure 3-8.- The ERTS-1 reflectance signatures of coastal features in MSS bands 5 and 
7 after corrections were made. (Data taken from all cluster analyses.) 
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3.3.4 Areal Measurement 

As a result of accurate discrimination of water bodies, an 
areal measurement technique was developed to estimate the area of 
lakes and ponds. To determine accuracy, areal measurements com- 
puted from the ERTS-1 data were compared with carefully measured 
areas computed from aircraft photography. An average error 
(i.e., the difference between area determined from ERTS-1 and 
aircraft data) in water body areal measurement was 2 square 
hectometers (5 acres) and uncorrelated with the size of the water 
body in the study area. The size of the lakes ranged between 
2 and 243 square hectometers (5 and 600 acres) and a variety of 
slopes . 

3.4 CONCLUSIONS 

The following conclusions can be stated. 

1. It was possible to use ERTS-1 data to map swamp features. 
Swamp areas were mapped to accuracies ranging from 33 to 96 per- 
cent. The lower accuracies were observed in swamp areas that had 
a significant forest cover. Such areas were classified as forest 
rather than swamp. This capability for wetland surveys could 
lead to better flyway management for waterfowl and a better 
understanding of the ecology of wetland systems. 

2. Computer-aided processing products were superior to 

those of conventional processing. The coastal features were 
surveyed by ERTS-1 in the following order of accuracy: water, 

salt marsh, grass vegetation, forest vegetation, marsh, nonvege- 
tated wetland, sand, and swamp. The ERTS-1 data are recommended 
for use as a basis for the specific application of detecting and 
mapping water bodies and coastal wetlands, for verifying circula- 
tion models for bays, and for measuring turbidity levels in 
estuarine water. 
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3. Ground surveys of turbidity in Galveston Bay strongly 
suggest that ERTS-1 data can be used to differentiate densities 
of undifferentiated particulate matter. In particular, channel 5 
(wavelength 600 to 700 nanometers) of the ERTS-1 photometer indi- 
cated a highly significant linear response to turbidity levels 
ranging between 20 and 120 p/m. Both supervised and nonsuper- 
vised classifications of ERTS-1 MSS data demonstrate relative 
differences in a radiance as a function of water turbidity in 
Trinity Bay. 

4. The ERTS-1 data could be used as a basis for surveying 
wetland areas to assist in management of tidelands. The areal 
extent of salt marsh areas was mapped to accuracies ranging from 
89 to 99 percent for three sets of ERTS-1 data (table 3-1). 

5. Long-term mapping of barrier island features (beach, 
urban, barrier flat vegetation zones, and bayside marshes) can 

be used to establish trends in barrier island evolution. Barrier 
island vegetation areas were mapped to accuracies ranging from 
87 to 100 percent for three sets of ERTS-1 data (table 3-II) . 

6. Water, low marsh, medium and high rrw :sh, nonvegetated 
wetlands, pine forest, deciduous forest, grassland, sand/concrete/ 
industrial-urban, and mixed vegetation and nonvegetated 

areas exhibited unique reflectance signatures for ERTS-1 data 
taken from August to November 1972, when atmospheric, Sun-angle, 
and instrument calibration corrections were made. 

7. Peak reflectance of water shifted from the green band to 
the red band for the highly turbid waters (particulate range of 
20 to 120 p/m) . 

8. Turbidity changes affected the IR reflectance as well as 
the visible reflectance of water for highly turbid water 

(>50 p/m) . 
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9. A phytoplankton bloom in Galveston Bay was discriminated 
best in ERTS-1 IR MSS band 6. 

10. Water coverage changes were detected from two sets of 
ERTS-1 data that were overlaid by digital processing means. 
Changes were related to flooding, tidal changes, and rice 
farming practices. 
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4.0 THE ERTS-1 FOREST ANALYSIS 

4 . 1 OBJECTIVES 

The forest investigation was performed < 1 } to determine the 
minimum size of forest features that could be detected in ERTS-1 
data under varying conditions, (2) to determine the suitability of 
the data for making forest classification maps by conventional and 
computer-aided methods, (3) to test the extension of feature clas- 
sifications from one area to another, and (4) to determine the 
accuracies of the classification results. 

4.2 ANALYTICAL APPROACH 

The Sam Houston National Forest in Texas (fig. 4-1) was 
selected as the test site because of the existing knowledge of the 
area and its accessibility to JSC for making field checks. The 
site comprises a primary study area surrounded by a secondary 
study area. The primary study area was used as a training area to 
develop methodology and techniques under closely controlled condi- 
tions. The secondary study area was used to verify those methods 
under conditions approximating an operational situation. 

Ideally, ground truth collected during the same time period 
as the acquisition of the ERTS data is necessary for a detailed 
evaluation of the data. The ground- truth map of the primary area 
(fig. 4-2) was prepared from available aerial photography col- 
lected during April, October, and November of 1972. The map was 
field checked on the ground and by helicopter. The ground- truth 
map of the secondary area was compiled after the area had been 
classified to ensure independent classifications. 
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Figure 4-2.- Ground-truth classification map of the primary area 

(NASA S-73-281 87) . 
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The hierarchy of forest land features established for this 
investigation (table 4-1) was based on expected ERTS-1 detection 
capabilities and on the needs of forest land management. The first 
three levels were established before ERTS-1 imagery had been re- 
ceived. On receipt of the data, it was apparent that the third 
level could be achieved and that a fourth level was a distinct 
possibility. Therefore, several Level XV classes were added to 
the hierarchy. This hierarchy of study features diverges from the 
USGS hierarchy (table 1-1) at Level II. The USGS scheme was 
established to fulfill the requirements of a general land use 
classification effort over the full gamut of possible landscapes. 
The scheme established for this investigation was restricted both 
to the requirements of forest management and to conditions that 
exist within the test site. 

The August 29, 1972, set of ERTS-1 MSS data was analyzed in 
this investigation. The data format evaluated included black- 
and-white film images from each of the four MSS bands, JSC color 
composites derived from digital tapes, GSFC color composites, and 
the digital data. 

The evaluation of ERTS-1 detection capabilities was performed 
with conventional image-interpretation techniques. The criteria 
established for the minimum size of features in the study area 
were as follows: 8 square hectometers (20 acres) for timber types, 

15 meters (50 feet) for width of streams and rights-of-way, and 
2 square hectometers (5 acres) for lakes and rural settlements. 
Thirty test targets (features) with known dimensions were estab- 
lished in the test site. Then, each of the four bands of black- 
and-white imagery, the JSC color composite, and the GSFC color 
composite were interpreted. The six image sources then were 
ranked, based on the number and size of test features detected 
(tables 4-II and 4-III) . 
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TABLE 4-1. - HIERARCHY DEVELOPED FOR THE FOREST INVESTIGATION 


Level II 


Level III 


Level IV 


Standing timber [Pine, established [Pine, established 

Hardwood, I Hardwood, 

established l established 

Cutover timber /Pine, cutover /Pine, cutover 

JPine, site Pine, site 

/ prepared prepared 

IPine, regenerated Jpine, site 

VHardwood, cutover | Prepared/ 

I vegetation 

Pine, regenerated 

Hardwood, cutover 


Streams and lakes! fstrearos and lakes 


Water 


Other land use 
classes 


Hardwood, 

established 

'Pine, cutover 

Pine, site 
prepared 

Pine, regenerated 
ardwood, cutover 


Impoundments 

rush 
Cultivated 
Grass 

Rights-of-way 
ural settlement 


[impoundments 

'Weeds 

Brush 

Cultivated 

Rights-of-way , 
dirt roads, or 
utility lines 

Grass 

Rights-of-way, 
paved roads 

Cut and bulldozed 

Bare soil 


Rural settlement 
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TABLE 4-II . - NUMBER OF TARGETS DETECTED 


Media 

Targets 

detected 

Percent 

detected 

JSC color composite 

21 of 30 

70 

Band 5 

20 of 30 

67 

Band 6 

16 of 30 

53 

GSFC color 
composite 

15 of 30 

50 

Band 7 

11 of 30 

36 

Band 4 

4 of 30 

13 


Classification of forest features involved both conventional 
image-interpretation and computer classification techniques. A 
black-and-white image of band 7, JSC and GSFC color composites, 
and multiband composites derived from the ACVP and the MCFV were 
classified using conventional interpretation techniques. Computer 
classifications included an unsupervised clustering algorithm 
(ISOCLS ) and a supervised classification procedure using a 
maximum- likelihood method (LARSYS). The basic classification pro- 
cedure involved establishing feature signatures in the primary 
study area, where 100-percent ground truth was used, then extend- 
ing these signatures to the secondary area, where ground truth 
was not compiled until the classification was completed. 

Accuracy of each classification was determined by comparing 
each classification map with ground truth, using two procedures. 

In the first procedure, each classification map was overlaid with 
a 0.64-centimeter (0.25 inch) grid that was also registered to the 
ground- truth map, and grid points of agreement between the two 
maps were counted. The total agreement point (correct classifica- 
tion points) divided by all points in the grid resulted in the 
? correct classification percentage for the classification. 
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TABLE 4-1 II.- 


MINIMUM SIZE OF FOREST FEATURES DETECTED 


Feature 

Smallest 

size 

2 

detected, hm 

Band or 
composite 
used 

Contrast 

with 

background 

Pine, established 

4 

5 

Good 



7 

Fair 



JSC 

Fair 

Pine, regenerated 

8 

5 

Poor 



JSC 

Fair 

Pine, site 

8 

5 

Fair 

prepared/vegetated 




Pine, cutover 

24 

5 

Good 



JSC 

Fair 

Hardwood, established 

12 

JSC 

Fair 

Brush 

4 

5 

Fair 



7 

Poor 

Grass 

6 

JSC 

Good 

Cultivated 

4 

5 

Good 



JSC 

Fair 

Streams and lakes 

2 

5 

Poor 



7 

Good 



JSC 

Poor 

Impoundments 

1 

7 

Fair 



JSC 

Poor 

Rights-of-way 

(a) 

5 

Fair 



JSC 

Fair 

Rural settlement 

40 

5 

Poor 


a 28 meters wide. 
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The second measure of accuracy for each classification map 
entailed computing the area of each class in the classification 
map and in the ground-truth map. The differences between the area 
totals enabled computation of an accuracy percentage. 

4.3 RESULTS 

The detection results are summarized as follows. Band 4, 
green (500 to 600 nanometers), was the poorest of the four bands. 
Band 5, red (600 to 700 nanometers) , was the best for most fea- 
tures but was poor for detecting hardwoods and water bodies. A 
pine stand as small as 4 square hectometers (10 acres) was detect- 
able when there was a good contrast with the surroundings. In 
band 6, near infrared (700 to 800 nanometers), hardwoods and water 
bodies were discernible, but pine stands diminished in detecta- 
bility because they faded into the background and mingled with all 
classes except hardwood and water. Band 7, near infrared (800 to 
1100 nanometers), had detection characteristics similar to band 6, 
with hardwoods and water bodies displayed even more prominently 
and other features displayed less prominently. The smallest fea- 
tures distinguished in band 7 were 1-square-hectometer (2.5 acre) 
ponds. These ponds were not distinguishable in the other bands or 
color composites. The composite image made at JSC produced the 
best detection results overall, and the GSFC composite rated 
between bands 6 and 7 (tables 4-II and 4-III) . 

4.3.1 Conventional Image Interpretation 

Only seven classes could be distinguished in band 7. The 
classification accuracy of the seven classes was 58 percent 
(table 4-IV). Ten classes of forest features were distinguished 
using the ACVP. The resulting map (fig. 4-3) was analyzed, and 
the classification accuracy by the point-count method was 55 per- 
cent and by the area-computation method was 64 percent. 
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TABLE 4-IV . - SUMMARY OF CLASSIFICATION ACCURACIES 


Technique 

Number 

of 

classes 

Classification 

a 

accuracy, 

percent 

Accuracy 
by grid- 
point 
count, 
percent 

Accuracy 
by area 
computation , 
percent 

Primary area 

Single band 

7 

50 

66 

58 

Optical multiband 

10 

72 

55 

64 

MCFV 

7 

50 

56 

82 

JSC color composite 

10 

72 

— 

67 

Clustering 

14 

100 

66 

81 

Maximum likelihood 

12 

86 

67 

73 

Secondary area 

Single band 

1 1 

79 

63 

83 

Optical multiband 

10 

72 

63 

84 

MCFV 


50 

56 

60 

JSC color composite 


72 

— 

78 

Clustering 

| 

100 

60 

75 


a Classif ication accuracies equal number of features identified 
divided by number of possible features (14) . 
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Figure 4-3.- Optical multiband classification map of the primarv 

area (NASA S-73-28195) . 



Seven forest features (fig. 4-4) were classified using the 
MCFV. Accuracy by the point- count method was 56 percent and by 
the area- computation method was 82 percent. 

Ten classes of features were distinguished by using the JSC 
color composite. These classes were annotated directly onto the 
DAS film recording (fig. 4-5), for which distortions had not been 
corrected. A skew distortion was introduced by the lag in each 
scan line, and a scale distortion was caused by different horizon- 
tal and vertical scales. Since no geometric correction was made, 
only the area method of accuracy computation could be used. 
Accuracy was 67 percent. 

4.3.2 Computer Classification 

Using nonsupervised techniques, 14 classes were recognized. 
The resulting printout map shown in figure 4—6 was converted to a 
cartographic format (fig. 4-7) at the scale of the ground-truth 
map. Accuracy by the point-count method was 66 percent. By the 
are a- computation method, accuracy was 81 percent. The most prev- 
alent error was the classification of regenerated pine as cutover 
hardwood. This error is understandable because the regeneration 
areas are cutover and the ground cover is a mixture of small pines 
and dense hardwood. 

Supervised classification techniques were also employed. 

First, a computer printout was produced by the LARSYS method. The 
printout was then changed to a cartographic format to match the 
ground-truth map. Accuracy by the point-count method was 67 per- 
cent. By the area-computation method, accuracy was 73 percent. A 
tape of the LARSYS classification was then processed on the DAS 
and recorded on film. This product is shown without geometric 
corrections in figure 4-8. 
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Figure 4-4.- The HCFV multiband classification map of the prima 

■rea (NASA S-73-28188) . 
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Figure 4-5.- The DAS multiband classification map of the primary 

area (NASA S-73-1 1 24) . 
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Figure 4-6.- Computer printout of cluster classification of the 
primary area (NASA S-7 j -2821 4) . 
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Approximate scale 1:63,360 
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Figure 4-7.- Cluster classification map of the primary area 

(NASA 3-7 3-281 89) . 
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The techniques previously described (except for the LARSY5 
classification) were applied to the secondary area using primary- 
area, spectral- signature data.. The accuracies achieved were com- 
parable to those in the primary area (table 4- XV) . The accuracy 
results in table 4- IV tend to support the feasibility and validity 
of extending signatures under the conditions of this investigation. 


Table 4-IV is a summary of the classification accuracies 
achieved in this investigation. Comparison of the results of the 
various classification procedures is difficult because not all 
parameters could be held constant. Clustering and maximum likeli- 
hood ranked higher than the conventional interpretation techniques 
in terms of the accuracy measures performed. 

The forest under investigation did not lend itself well to 
precise classification. Timber classes occurred in mixtures of 
species, age, vigor, and size, blending gradually from one to 
another. Even under these conditions, the are a- computation accu- 
racies achieved were quite high. Table 4-V contains the class 
(feature) accuracies achieved in the primary study area. Both in 
terms of the number of features classified and individual feature 
accuracies , the computer-aided classification techniques were 
better than the conventional image- interpretation techniques. 


Other items of interest to foresters , but incidental to the 
investigation, were found in the ERTS-1 data. For example, alight 
ground fire, a prescribed burn of approximately 40 square hecto- 
meters (100 acres) to clear brush, was discovered. A JSC color 
composite of this coverage (fig. 4-9) shows the burn as a rusty 
brown area. The precision with which the effects of this light 
fire were registered indicates ■chat ERTS— 1 data could be used in 
fire-damage assessment to map the perimeter of large burns. 


TABLE 4-V. — ACCURACY OF FOREST FEATURE CLASSIFICATIONS BY THE ARE A- COMPUTATION METHOD IN THE 

PRIMARY STUDY AREA 



Accuracy, percent 


Single Optical MCFV JSC color Clustering Clustering Maximum 
band multiband multiband composite (small area) (large area) likelihood 
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Pine, cutover, and hardwood, cutover, were combined into a single class 
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Figure 4-9.- The DAS enhancement of light ground fire in the Sam 
Houston National Forest (NASA S-73-25467) . 
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During the investigation, a pine bark beetle infestation 
occurred vrithin the test area. Damage was detected indirectly in 
,1 one case? pines covering an area of approximately 2 square hecto- 

meters were killed in a stand classified as hardwood, although 
i nearly 50 percent of the stand was pine. The L&RSYS program 

i classified this portion of the stand as cutover hardwood. Because 

LARSYS is a supervised program, a beetle-infested training field 
would, have to be included to enable investigation of such a 
' j classification. 

4 . 4 CONCLUSIONS 

The following conclusions can be stated. 

1. The ERTS-1 data can probably best be used in forestry 

applications if the data are used for extensive surveys in which 
broad, generalized classes are needed, rather than for intensive 

1 surveys in which detailed stand conditions must be portrayed. 

2. The clustering and maximum- likelihood methods of classi- 
fication are both efficient, and differences in accuracies achieved 
are not significant enough to rank one above the other. However, 
computer methods as a group can probably be ranked as superior to 
conventional methods for forest classification. 

3 . in addition to the possibilities of producing timber 
inventory maps, the ERTS-1 data have the advantage of providing 
sequential coverage and thus providing an increased capability to 
classify certain categories such as evergreen, deciduous, and 
mixed. 
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5.0 THE ERTS-1 RANGE ANALYSIS 
5.1 OBJECTIVES 

The primary objective of the ERTS- 1 Range Analysis Team 
investigation was to determine the usefulness of ERTS-1 MSS data 
in classifying types of vegetation for mapping range and related 
grazing lands . Assessment of the capability to use ERTS-1 MSS 
data to detect, identify, locate, and determine the area of range- 
land vegetal types was therefore required, i. 

In USGS Circular 671, Anderson and. others define the land use 
category "rangeland" as being composed essentially of grass and 
shrub vegetation types. Current management practices define 
rangeland more broadly, as consisting of both nonwoodland and 
woodland vegetation types. In this study, the latter definition 
was adopted and many species of trees, brush, grasses, Shrubs,, and 
forhs associated with the two broad categories were considered. 

The rangeland classification hierarchy adopted by the ERTS-1 Range 
Analysis Team to meet the primary investigation objective is shown 
in table 5-1. 

Specific objectives in the classification of rangeland sub- 
groups reflected the broad definitions of rangeland described 
previously and included the following. 

1. To classify the Level II rangeland classes of woodland 
and nonwoodland, as opposed to nonrangeland such as cropland, 
water, and urban areas 

2. To classify Level III classes 

3. To consider an assessment of 

a. The accuracy of vegetation type identification 

b. The precision of areal measurements from ERTS-1 data 

Co The use of temporal signatures in classification 


TABLE 5-1.- RANGELAND CLASSIFICATION HIERARCHY 


Level I; Level II Level II I a 

("Post oak stands 

r 

Woodland * Mesquite stands 

^Bottom land hardwood 

Rangeland ■< Wetland, marshhay cordgrass 

Gulf cordgrass 

Nonwoodland H Smutgrass 

Burned gulf cordgrass 
Cultivated bermuda- fields 

- - (Abandoned cropland 

Deep seawater 
f Water * Coastal water 

Nonrangeland i Inland Water 

Urban areas ) Residential and the like 

I ^ Cropland j | Rice fiel ds and the like 

a Al though not necessarily representative of the rangeland 
features m the United States, Level -III features in this table 
were chosen because they were dominant in the study , sites of 
this investigation. Since nonrangeland was not of interest in 
tnis study. Level III categories were not defined in detail. 


Level II 


Rangeland 


Nonrangeland 


Woodland 


Nonwoodland 


Water 


Urban areas 



5.2 ANALYTICAL APPROACH 


Two study sites in the Houston Area Test Site were selected 
for the rangeland investigation. These sites were designated the 
Snook site and the San Bernard site. The Snook site occupies the 
boundary between the Texas black-land prairie region and the Texas 
claypan savannah. Features of interest in the Snook site are 
native stands of post oak, bottom land hardwood, mesquite, native 
grassland, and planted pastures of bermuda grasses. The San 
Bernard site is typical of the gulf coast marsh region along the 
Texas and Louisiana coasts. Two distinct vegetation zones existed 
in this site, a marshhay and smooth cordgrass zone (the salt marsh 
zone) and a gulf cordgrass zone (the slightly higher zone). Both 
sites offered the opportunity for Level II classification, whereas 
the San Bernard site provided a potential for Level III classifi- 
cation of vegetation types. 


Conventional image-interpretation techniques and computer- 
aided data-processing techniques were used for analyzing MSS 
imagery and digital tape data. Color enhancements were obtained 
by conventional film data processing. Three photointerpreters 
who were not familiar with the study areas were asked to delineate 
the features of interest in the enhancements to obtain unbiased 
interpretations . The delineated features then were planimetered 
three times and the areal measurements averaged. Ground-truth 


information from low— altitude-aircraft photography was used to 


check the measurements obtained. A regression analysis also 
made to investigate the significance of the results.. 


was 


5 . 3 RESULTS 


A temporal color enhancement of the Snook site, created on an 
MCFV , is shown in figure 5-1. Figure 5-2 is a classification map 
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Figure 5-2.- Classification map of the San Bernard si 
supervised computer classification techniques (ERTS 
1972, digital data) (NASA S-73-28326) . 
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of the San Bernard site made using supervised computer classifica- 
tion techniques. Figure 5-3 is a classification- like clustering 
map of the Snook site made using nonsupervised computer classifi- 
cation techniques. Figure 5-4_.is a representative bar-chart 
illustration of the computer classification accuracies for the 
Snook single-date and temporal data. The team also performed a 
linear regression analysis of areal measurements obtained by plan! - 
metering color enhancements of ERTS-1 film imagery. 


5.3.1 Conventional Data Processing 

The following was accomplished or determined using conven- 
tional data processing. 

T. An examination of the color enhancements created on the 
film viewers indicated that Level IX classification of both sites 
appeared satisfactory. The sites were separated into wooded 
rangeland, nonwooded open rangeland, cropland, water, and urban 
areas. 

2. The Level III features of wet lowland marshhay were 
separated from the drier upland zone of gulf cordgrass in the San 
Bernard site. 

3. Regression analyses were performed to correlate the areas 
interpreted from color enhancements with ground truth. A linear 
model was developed that would predict the true areal measurement 
using a measurement obtained from a color enhancement created on 
available film viewers. However, the errors of the areal measure- 
ments were found to be as high as 214 square hectometers 

(530 acres) when measuring areas as large as 9100 square hecto- 
meters (22 500 acres) . The inordinately large error figures were 
attributed to the fuzziness of boundaries between the extremely 
complex vegetation zones , and to the less than ideal quality of 
the color enhancements created. 



ES Vigorously growing cropland 
I i Water 
□ Clouds 

H Rangeland, with woodland 
I i Rangeland, open rangeland, and grassland 


gure 5-3.- Classification-like clustering map of the Snook site 
made using nonsupervised computer classification techniques 
(ERTS— 1 August 30, 1972, digital data) (NASA S— 73— 28028) . 
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Figure 5-4.- Snook data overall classification accuracy. Snook 1 
and Snook 2 ERTS-1 data were taken August 30 and November 28, 
1972, respectively. The Snook 1/2 category is the temporal 
composition of Snook 1 and Snook 2 data. 
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4. Some features (e.g. , cropland) were enhanced and accented 
from the background more sharply In temporal enhancements than in 
single-date enhancements. 

5.3.2 Computer-Aided Data Processing 

The following was determined using computer-aided data 
processing. 

1. Classification accuracies of training fields were high. 
When the fields classified at Level III were aggregated into 
Level II categories, the Level II classification accuracies 
appeared even higher, as high as the mid- 90 -per cent range 

(fig. 5-4) . 

2. Temporal analyses gave even better classification accu- 
racies for the training data. Furthermore, confusion between 
training classes was reduced. 

3. Because of the lack of intensive ground-truth data and a 
deficiency of spot checks, many test fields were not statistically 
representative. Consequently, classification accuracies of test 
fields were not as high as those of training fields. Level II 
classification accuracies appeared higher than those for Level III. 

4. In general. Level II classification was satisfactory for 
both sites (i.e. , the separation into wooded rangeland, nonwooded 
rangeland, cropland, and water) . In the San Bernard site, the 
soil moisture content of marshhay, cordgrass, and gulf cordgrass 
was sufficiently different to permit separation between the two 
classes . 

5. Clustering results were very useful because no a priori 
information was required to identify spectrally homogeneous areas 
in the data. Furthermore, this nonsupervised technique allowed the 
grouping of features for which ground truth was not available into 
spectrally unique groups. Further analysis could, in many cases, 
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result in assigning significance to these classification-type 
cluster maps. An attempt to use a sampling technique to reduce 
the computation time of the clustering process also was made. 

6. Clustering results appeared to be very satisfactory for 
Level XI classification and, to some extent, for Level ILL classi- 
fication. Because of a lack of extensive ground truth, a quanti- 
tative analysis of the classification accuracies was not made. 

5.4. CONCLUSIONS 

The following conclusions about the rangeland study resulted 
from an analysis of ERTS-1 MSS data using conventional image- 
interpretation and computer-aided processing methods. 

1. The ERTS-1 Range Analysis Team demonstrated the useful- 
ness of ERTS-1 MSS data for satisfactory classification of vegeta- 
tion types into broad Level II categories. Rangeland comprising 
woodland and nonwoodland was distinguished from nohrangeland com- 
prising water, urban areas, and cropland. 

2. Finer Level III classification also was possible but 
depended on the existence of certain characteristics of features 
such as soil moisture content that are not yet fully understood. 
The demonstrated capability to separate the wet lowland zone from 
the drier upland zone (e.g. , in the San Bernard site) should be of 
significant value to users interested in mapping wetland features 
in coastal areas. 

3. Temporal analysis, particularly using computer-aided 
techniques xn processing digital data, resulted in enhanced dis- 
crimination between nearly all vegetation types encountered in 
this investigation. This conclusion is in accord with the predic- 
tion from mathematical theories. 
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6.0 THE ERTS-1 LAND USE ANALYSIS 

6.1 OBJECTIVES 

The general objective of the ERTS-1 land use investigation 
was to evaluate how well data from the ERTS-1 MSS could be used 
to detect, identify, and delineate land use features within the 
HATS, an 18-county area around Houston established previously as 
a land use test area. A more specific objective was to determine 
whether the land use classification scheme proposed in USGS Cir- 
cular 671 could be used as the basis for delineating land use 
by conventional image interpretation and computer-aided classi- 
fication of ERTS-1 data. 

6.2 ANALYTICAL APPROACH 

An analysis of the entire 41 000-square-kilometer 
(16 000-square-statute-mile) HATS area was not feasible with the 
available man-hours and computer time allotted to this investiga- 
tion. Consequently, a 4660-square-kilometer (1800-square-statute- 
mile) study area was selected to correspond to the data on one 
CCT . The study area is oriented generally north-south and repre- 
sents one-fourth of a scene of ERTS imagery (fig. 6-1). 

An attempt was made to delineate Level I land use categories 
as urban and built-up land, agricultural land, rangeland, forest 
land, nonforested wetland, water, and barren land by conventional 
image-interpretation techniques. Black-and-white images of the 
study area from ERTS-1 MSS bands 5 and 7 (October 4, 1972, pass) 
were enlarged to a scale of approximately 1:250,000, and deline- 
ations of land use were recorded on transparent overlays of these 
enlargements. A comparative study was conducted by using similar 
interpretative techniques to delineate land use categories on 






6-3 


enlargements made from first-generation color composites obtained 
directly from a DAS film recorder at JSC. These enlargements were 
also at a scale of approximately 1:250,000 and were simulated 
color-IR composites of bands 4, 5, and 7. The same techniques 
were used in delineating some Level II land use categories. 

Two basic computer-aided classification techniques (super- 
vised and nonsupervised) were used in classifying the study area 
into land use categories. The ISOCLS, a nonsupervised clustering 
algorithm, was used to group every sixth picture element (oixel) 
from every sixth scan line into clusters of pixels having similar 
spectral characteristics. This reduction in the number of data 
points (40 000 pixels) was necessary because the capacity of the 
computer was not sufficient to process the total number of data 
points (1.3 million oixels) covering the entire study area. This 
3-percent systematically alined sample of data points, uniformly 
distributed over the entire study area, was grouped into spectral 
clusters, each representing a portion of the full range of spec- 
tral variations found in the study area. The input parameters to 
the cluster program could be adjusted to provide more clusters or 
fewer clusters. However, after considering the amount of detail 
needed for the proposed land use hierarchy and the estimated com- 
puter time required, it was reasoned that input parameters provid- 
ing 13 clusters would be a reasonable compromise. Graymap 
printouts depicting the spatial distribution of the pixels grouped 
into each cluster were generated on the computer. Each cluster 
was identified and assigned to a specific land use category by 
correlating the cluster delineations on the graymaps with exist- 
ing land use maps and aerial photographs, and by analyzing perti- 
nent cluster statistics that had been plotted on graphs. After 
grouping the clusters into the desired land use categories, a 
color-coded cluster map in the form of a color transparency was 
produced on the JSC DAS film recorder. 


I 
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Once the clusters had been grouped satisfactorily into the 
Level I land use categories, the means and covariance matrix 
statistics from the cluster analyses were substituted for training 
field statistics as inputs in the LARSYS II supervised classifica- 
tion approach. The use of cluster statistics in lieu of training 
field statistics eliminated some of the difficulties that would 
have been encountered in selecting representative training fields 
for such a large study area for which intensive ground truth or 
large-scale aerial photography was not readily available for anal- 
ysis. Because of the relative spectral complexity of much of the 
study area landscape, it was deemed desirable to be able to clas- 
sify every pixel (instead of every sixth pixel) within the entire 
study area. To do so, however, it was necessary to divide the 
entire area into north-south linear strips, with the number of 
data prints in each strip not to exceed the storage capacity of 
the computer memory drum. Level I land use classification maps of 
each strip were subsequently mosaicked to form a classification 
map of the entire study area. Experience gained in delineating 
Level I land use categories by both supervised and nonsupervised 
classification techniques indicated that a potential existed for 
dividing the urban and built-up category into some Level II cate- 
gories. Because some urban features (vegetated residential areas) 
have spectral characteristics similar to some nonurban features, 
such as forest and agricultural areas, it was necessary to reas- 
sign the 13 original clusters to land use categories that would 
represent three Level II urban and built-up categories (residen- 
tial, commercial/industrial/transportation, and open) when a 
Level II classification was made of only the urbanized portion of 
the study area. 

The accuracy of the three classification approaches was 
assessed by measuring the agreement between the classified data 
and the base reference data established for the accuracy analysis. 
Five accuracy test sites, ranging in size from 21 to 104 square 
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kilometers (8 to 40 square statute miles), were established m the 
study area. Base reference data were established by visually 
classifying land use in each accuracy test site from high- 
tude, infrared-Ektachrome photography acquired April 22, 1972. 

Each site was divided into 2 . 6 -square-kilometer quadrats, and the 
percent occurrence of each class in each quadrat was measured 
using a dot-sampling technique. The same procedure was performed 
on each class for each classification product except for the 
computer-aided classification maps, for which pixels in each 
class were counted and converted to percent occurrence. The per- 
cent agreement (class-by-class comparison of accuracy) between 
classification products and base reference data was calculated on 
the basis of the percent occurrence. 


6 . 3 RESULTS 


The following results were obtained. 


1 . a visual comparison of all the classification results 
(figs. 6-2 to 6-7) shows a strong correlation in the areal pat- 
terns of land use between the three analysis approaches used in 
the investigation. However, there is a significant difference m 
detail. Because of the relatively small scale (1:250,000) of the 
manually interpreted imagery, many of the smaller features were 
difficult to portray. The result is a pattern of relatively 
homogeneous tracts of land use classes (figs. 6-2 to 6-3). 


2. The computer-aided classification maps display a finer 
texture in the land use patterns (figs. 6-5 to 6-7). This finer 
precision is a result of the capability of the computer to clas- 
sify each pixel (approximately 0.45 square hectometer). 

3. The image interpreter can compensate for his inability to 
resolve fine details with the ability to resolve spatial patterns 
and relationships in the land use features. This capability was 
particularly prominent in the urban areas, where many linear 
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Figure 6-2.- Level I land use classification of study area (ERTS-1 
bands 5 and 7 from October 4, 1972) (NASA S-73-31699). 
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Figure 6-3.- Level I land use manual interpretation of study area 
(ERTS-1 bands 5 and 7 from October 4, 1972) (NASA S-73-25517). 
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Figure 6-4.- Level I 
(ERTS-1 bands 4, 5 


and II classification of land use study area 
and 7 from August 29, 1972) (NASA S-73-28191) 
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Figure 6-5.- Land use cluster nap developed fron a 3-percent data 

sample (NASA S-73-28073) . 






Urban/built-up 
CZI Agriculture/rangeland 
IZZ] Forest 
Water 

CU Wetland 


Computer classification of Level I land use 
(NASA S-73-2853) . 







Jczj Residential 

HU) Commercial, industrial, and roads 
CZ3 Open and other 
CT3 Agriculture/range 
era Forest 


Wetland 


Computer classification of land 
(urban) (NASA S-73-2866) 


ORIGINAL PAtife fe 


6-12 


features (secondary roads for instance) could be visually distin- 
guished by conventional image interpretation even though the 
widths of the features were well below the spatial resolution 
threshold of the scanner. 

4. Relatively high classification accuracies for Level I 
land use categories (table 6-1) were achieved by conventional 
image interpretation and computer-aided classification techniques, 
with the exception of the urban and built-up category when it was 
derived from computer classification of the entire study area. 

When only the preselected urban area was classified in Level II 
categories, considerably better computer classification accuracies 
were attained (table 6-II) . This apparent discrepancy in accu- 
racies was probably due to the spectral heterogeneity of the urban 
scene which also contained vegetated agriculture-rangeland 
features . 


6.4 CONCLUSIONS 

The following conclusions can be stated. 

1 . It was concluded from this investigation that general 
land use categories as suggested for Level I and some Level II 
categories in USGS Circular 671 could be obtained over relatively 
large areas from ERTS-1 MSS data by conventional image interpre- 
tation and computer-aided classification techniques. 

2. In the computer-aided processing, a small (3 percent) 
sample of the available digital data was sufficient to identify 
the general land use categories throughout the entire study area. 
This conclusion indicates that even larger geographic areas could 
be similarly classified without exceeding nominal computer 
capacities . 

3. Where greater classification accuracies or more detailed 
land use categorizations of larger areas are desired, it may be- 
come necessary to define categories of land use by geographic 
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TABLE 6-1.- LEVEL I AND II LAND USE PRODUCT AGREEMENT WITH BASE DATA 

[Based on percent occurrence] 


Land use class 


Base quadrats 

Sample Number 

test (100 Class occurrence 

site point 

counts 

each) Count Percent 


Conventional image interpretation 


Forest 


Cumulative total 
Agriculture/ range land 

Cumulative total 
Water 

Level I urban 

Level II urban 

Residential 

Commercial/ 

industrial/ 

transportation 

Open and other c 


325 54.2 


a ' b 39 a ' b 3724 a,b 95.5 


Black-and-white 
imagery class 
occurrence 



Count 

Percent 

Count 

2394 

72.5 

2800 

3473 

91.4 

3455 

293 

19.5 

265 

267 

66.8 

286 

6427 

71 .4 

6806 

552 

32.5 

■8 

100 

14.3 


2981 

85.2 


3633 

61.6 

3734 

259 

43.2 

247 

b 1 100 

b ioo.o 

3855 


Color composite 
imagery class 
occurrence 


Percent 




Computer 
classification 
class occurrence 

Count 

Percent 

2383 

72.2 

3340 

87.9 

240 

16.0 

192 

48.0 

6155 

68.4 

719 

42.3 

214 

30.5 

2825 

80.7 

3758 

63.7 

296 

49.3 

1642 

42.1 

1514 

n 

B2 

m 

124 

10.3 


a August 29, 1972, ERTS-1 data used for all analyses except those using October 4, 1972, data. 

b 0ctober 4, 1972, ERTS-1 data base covered only part of test site in black-and-white imagery. For base 

quadrats, number = 11; class occurrence count = 1035; and class occurrence percent = 94.1. 

c Category not delineated by conventional image interpretation. 
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TABLE 6-II . - LEVEL II URBAN LAND USE PRODUCT AGREEMENT WITH BASE DATA 


[Based on percent occurrence] 


Land use 
category 

Classification product 

Sample 

size 

JSC color 
composite, 
percent 

Sample size, 
percent 

Computer 

aided, 

percent 

Base 

1 Counted 

1 _ 

Base j 

Counted 

Residential 

1293 

mm 

64.7 

8£1M 


83.0 

Commercial/ 

1 37 


43.8 

■ 


60.2 

industrial/ i 




I 



transportation ; 




I 



Open and other 

— 


— 

1 

10.3 

61 .7 
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region, perform sampling within each region, and classify the 
entire large area into the desired land use categories using 
computer-aided techniques. 
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7 . 0 THE ERTS-1 URBAN LAND USE ANALYSIS 

7.1 OBJECTIVE 

The objective of the ERTS-1 urban land use investigation was 
to evaluate how well data from the ERTS-1 MSS could be used to 
detect, identify, and delineate urban features within the Houston, 
Texas, metropolitan area. 


7.2 ANALYTICAL APPROACH 


Because of the complexity of urban landscapes and the esti- 
mated amount of computer time required, the scope of the inves- 
tigation was restricted to four separate residential study areas 
and to the major transportation network around Houston. During the 
computer-aided classification phase of the investigation, it was 
found advisable to expand the study area to include more varieties 
of urban land use. A larger, contiguous area embracing approx- 
imately the eastern two-thirds of the Houston metropolitan 
area and encompassing the four residential study areas was selected 
for the final investigation. 

The set of ERTS-1 MSS data covering the subject area 
obtained August 29, 1972, was used in this investigation. Appro- 
priate frames of aerial photography from April and October 1972 
were used as collateral ground-truth data. The preliminary urban 
land use classification scheme used in this investigation was a 
modification of the land use scheme proposed in USGS Circular 671 
{table 1-1). 

Three basic analytical approaches were used in attempting to 
meet the objectives of this investigation: conventional image 

interpretation, correlation of computer graymaps and aerial 
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photography, and computer-aided classification. Several differ- 
ent analytical techniques were devised to evaluate the usefulness 
of these approaches. 

Conventional image-interpretation techniques were used to 
analyze ERTS-1 enlarged black-and-white imagery and color com- 
posite imagery generated from the digital data. Supervised and 
nonsupervised computer-aided classification techniques were used 
with ERTS-1 computer-compatible tapes. 

Two data-processing procedures were used in the supervised 
computer-aided techniques. In one procedure, the Earth Resources 
Interactive Processing System developed at JSC was used in con- 
junction with the LARSYS II classifier on the IBM 360 series com- 
puter. In the second procedure, the Module Training Field Option 
(MTFO) was used on the Univac 1108 series computer. A LARSYS II 
classifier also was used in this procedure, but the investigator 
could modify the statistical inputs. 

The nonsupervised computer-aided classification investiga- 
tions involved the use of ISOCLS clustering to examine the spec- 
tral composition of picture elements of clusters within a complex 
urban scene, and to determine the accuracy to which the various 
clusters could be correctly grouped to identify specific urban 
land use categories. The accuracy analyses of the various 
computer-aided classification techniques were performed by corre- 
lating samples of individual picture elements to specific 
geographic areas as identified and delineated on aerial photog- 
raphy by a picture-element-correlation grid technique developed 
during this investigation. 

A separate attempt was made to delineate the major trans- 
portation routes in the Houston metropolitan area as a major 
urban land use category . Conventional image interpretation 
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techniques were used in delineating the interstate freeway system 
and some of the major highway routes and streets from enlarged 
ERTS-1, band 5, black-and-white imagery (scale approximately 
1:250/000). Both the LARSYS and the ISOCLS computer-aided classi- 
fication techniques were also used to produce maps of the major 
transportation routes. Two approaches were used in classifying 
with the ISOCLS technique. In one approach, clustering based on 
ERTS-1 bands 5 (red) and 7 (infrared) was used. In the second 
approach, clustering based on ERTS-1 bands 4 (green) , 5 (red) , 

6 (infrared), and 7 (infrared) was used. 

7.3 RESULTS 

In the conventional image-interpretation investigations, the 
importance of using spatial pattern recognition clues for inter- 
preting the extremely small-scale ERTS-1 imagery was emphasized. 
The following results were obtained. 

1 . The spatial resolution of the multispectral scanner was 
sufficient for recognizing only gross geographic patterns, rather 
than any detailed textures that could provide clues to the iden- 
tity of certain urban features. However, the limited spatial 
resolution and the extremely small scale of the imagery did com- 
bine to present a textural pattern characteristic of the highly 
built-up areas of much of metropolitan Houston, where the wide 
streets and rows of bright rooftops gave a distinctive cross- 
hatched texture to the imagery. Linear patterns indicating the 
major highways and streets were also readily recognizable, 
despite the fact that the widths of these features were well 
below the spatial resolution capability of the scanner. It was, 
therefore, possible to differentiate visually certain urban 
features from the surrounding nonurban landscapes by interpreting 
the extensive linear and crosshatched patterns as spatial surro- 
gates of urban features. 
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2. Dif f iculties were encountered in manually delineating 
some recognizable urban features simply because of the extremely 
small physical dimensions, even when the imagery was photograph- 
ically enlarged to its limit. 

3. As the initial results from the computer classification 
experiments were being compared to ground-truth data, it became 
apparent that considerably more detailed urban land use informa- 
tion was surfacing than had originally been anticipated. 

4. The residential land use category was actually one of the 
most complex and least consistent categories to be delineated by 
spectral classifications. The residential category was a spec- 
trally heterogeneous intermixture of small vegetated and non- 
vegetated surfaces. For this reason, the study areas were 
expanded to include a greater variety of urban land use categories 
so that a broader statistical base would be available and greater 
spectral contrasts would be available, in addition to those 
adjacent to the residential areas. 

5. The accuracy to which the original selection of residen- 
tial study areas could be differentiated from the surrounding land 
use categories by means of computerized classification techniques 
is shown in tables 7-1 and 7-II. Table 7-1 contains a comparison 
of three computer classification techniques used in classifying 
various urban land uses associated with selected residential 
study areas. The ERIPS and MTFO techniques were the two 
approaches used in the supervised classification programs. In the 
ISOCLS technique, nonsupervised classification programs were used. 
The accuracy of these computerized techniques for differentiating 
the selected residential areas from adjacent land use areas was 
calculated by determining the proportion of the arbitrary bound- 
ary that could actually be differentiated. 

6. Another measure of accuracy was obtained by comparing 
the acreage contained in each residential study area as measured 
on aerial photography with the acreage determined from the ISOCLS 
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TABLE 7-1.- ACCURACY OF DIFFERENTIATION OF LAND USE BOUNDARIES BY 
COMPUTER-AIDED CLASS T FICATION TECHNIQUES 


Residential 
study area 


Percent of boundary length 
correctly differentiated — 

Land use category 

By ISOCLS 
techniques 

By ERIPS 
techniques 

By LARSYS 
techniques 

Cloverleaf 

Major transportation 

0 

58 

28 

(area L) 

Water (stream) 

100 

(a) 

(a) 


Forest 

56 

50 

50 


Open fields 

32 

(a) 

(a) 

Pasadena (area M) 

Commercial/ industrial 

(a) 

(a) 

68 

Garden Villas 

Major transportation 

59 

(a) 

(a) 

(area N) 

Water (channel) 

26 

(a) 

(a) 


Open fields 

48 

(a) 

(a) 


Other residential 

(a) 

34 

30 

Center City 

Commercial/ industrial 

34 


27 

(area P) 

Major transportation 

57 


(a) 


Institutional 

78 

(a) 

(a) 


Water (channel) 

24 

(a) 

(a) 


Open fields 

22 

(a) 

(a) 


Other residential 

(a) 

0 

18 


a No comparative analysis was made. 


















7-6 


TABijE 1 - 11 .- ACCURACY OF AREAL EXTENT OF LAND USE CATEGORIES 
DELINEATED BY COMPUTER-AIDED CLASSIFICATION TECHNIQUE ISOCLS 


Residential study area 

Ratio of area acreage C/B a 

Cloverleaf (area L) 

1.12 

Pasadena (area M) 

(b) 

Garden Villas (area N) 

1 .16 

Center City (area P) 

.91 


B _ acreage of area as delineated on an aerial photograph; 

C = acreage of area as delineated on computer classification map. 

No comparative analysis was made. 




1 


7-7 

output maps. Table 7-II includes the ratio of areas calculated 
by the ISOCLS technique. A ratio greater than 1.0 indicates that 
the computer classification technique overclassified the size of 
the area, whereas a ratio less than 1.0 shows that the area 
classified by the computer was smaller than the actual size of 
the area. 

7. Determining classification accuracies of areas encom- 
passed by boundaries arbitrarily selected from panchromatic aerial 
photographs may have imposed unreasonably stringent requirements 
on the classification techniques. For practical reasons, a com- 
promise was necessary between the amount of digital data that 
could be processed within a reasonable turnaround time and the 
number of residential study areas having portions of boundaries 
that may not be entirely distinguishable, even by ground obser- 
vation. To determine the accuracy of the land use map generated 
by using nonsupervised classification, a misclassif ication error 
matrix was generated for the categories developed from the 
32 clusters (table 7-III) • The misclassif ication error matrix was 
established by comparing the nonsupervised classification with a 
base developed from 1 : 1 20 , 000-scale high-altitude-aircraft pho- 
tography (ground truth) of the same area. Table 7-III includes 
the overall accuracies obtained when the residential land use 
category was not limited to arbitrary boundaries but was inte- 
grated with other land use categories over a major portion of the 
entire Houston metropolitan area. 

8. Figures 7-1 to 7-3 provide an overview of the complex- 
ities of the urban land use categories delineated by the ISOCLS 
classification techniques and the two LARSYS classification 
techniques. A cursory comparison of the three figures indicates 
an impressive overall similarity of urban land use patterns. A 
more detailed examination of figures 7-1 to 7-3 discloses some 
differences in the areal extent of certain categories. Most of 
these differences can be attributed to the manner in which a 


TABLE 7-III.- MISCLASSIFICATION ERROR MATRIX 


[Percent agreement with base] 


Category 

Commercial/ 

industrial/ 

transportation 

Residential 

Urban 

Vegetation 

(woody) 

Vegetation 

(nonwoody) 

Water 

Commercial/ industrial/ 

94.2 

5.5 

— 

— 

— 

0.3 

transportation 

Residential 

2.6 

66.8 

23.0 

4.5 

3.2 

— 

Mixed urban 

1.0 

20.8 

51.1 

3.8 

23.5 

— 

Vegetation (woody) 

— 

.7 

.2 

95.1 

4.0 

— 

Vegetation (nonwoody) 

1.1 

12.1 

25.7 

4.8 

56.2 

— 

Water 

3.9 

3.0 

1.9 

2.2 

1.5 

87.7 
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Figure 7-1.- Urban nonsupervised computer-aided classification 
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Figure 7-2.- The ERIPS supervised computer-aided classifications 
delineating urban land use (NASA S-73-25486) . 
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Figure 7-3.- 


The MTFO supervised computer-aided classifications 
delineating urban land use. 
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specific spectral cluster was assigned to a particular land use 
category. In the supervised classification, certain differences, 
especially in the residential category, could be attributed to 
difficulties in finding specific training areas representative of 
that class of features. Ellington Air Force Base, located near 
the lower right corner of each figure, illustrates consistency 
among the three techniques. Note that the area to the left 
that contains barracks and other buildings is classified as resi- 
dential, whereas the hangars and runways on the right are classi- 
fy ed as conunercial/industrial/ transportation . Some vegetated 
areas are indicated between the runways. 

9. Only a qualitative study was made to evaluate the use- 
fulness of ERTS-1 MSS data for delineating the major transporta- 
tion routes within the metropolitan Houston area. Maps produced 
both by conventional image interpretation (fig. 7-4) and by 
cluster-map computer classification procedures (fig. 7-5) were 
compared with a base highway map of Houston (fig. 7-6). This 
comparison revealed that major freeways, highways, and many major 
Houston streets could be delineated from the ERTS-1 imagery and 
digital data. This result suggests that automated computer pro- 
cedures could be used with a very high degree of confidence in 
delineating a substantial portion of the major highways and 
street pattern., in a complex urban area. 

10. The investigations involving the correlation of com- 
puter graymaps with aerial photographs revealed that relatively 
unsophisticated conventional image-interpretation techniques 
could play an important role in exploiting the full capability of 
ERTS-1 data in urban land use analyses. The use of specially 
constructed reference grids demonstrated that the graphical posi- 
tion of a picture element could be correlated from a computer 
printout map to its precise location on an aerial photograph. 

This procedure made possible interpretation of the meaning of the 
classification anomalies encountered in the computer 
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Figure 7-4.- Major transportation routes in the Houston area. 
(Conventional image interpretation from an ERIPS supervised 
classification of ERTS-1 MSS data, August 29, 1972.) 


Figure 7-5.- Cluster map of urban residential area showing high 
reflectance areas from nonsupervised classification 
August 29, 1972. 







Base highway map of Houston. Heavy outline (lower 
right) encloses area shown in figure 7-5. 
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classification programs, and also aided in evaluating the accura- 
cies of the various outputs from these programs. 

11. The initial guidelines for this investigation stipulated 
the USGS proposed land use scheme as the goal for classifying 
urban land use. This conventional land use scheme was not appro- 
priate when the computer classification programs had ERTS-1 
spectral data as the only input. Table 7-IV is the classification 
hierarchy developed on the basis of spectral contrast levels that 
appeared to correlate most consistently with the capability of the 
various computer classification programs. Some shortcomings are 
recognized in trying to adapt this land use classification scheme 
to a user-oriented requirement, and it is offered here only as a 
suggested point for discussion when comparisons can be made with 
other land use classification schemes devised by ERTS-1 inves- 
tigators of urban areas in other geographic regions. 

TABLE 7-IV.- COMPUTER-AIDED CLASSIFICATION SCHEME 
FOR LEVEL I AND II URBAN LAND USE 


Level I 

Level II 

Nonvegetated 
land areas 

Commercial/industrial/ 

transportation 

All other 

Residential 

Vegetated 
land areas 

urban 

Mixed urban 

Nonwoody vegetation 


Woody vegetation 

Water 

Water 

Clouds 

Clouds 

Shadows 

Shadows 
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7.4 CONCLUSIONS 

The following conclusions can be stated. 

1 . The ERTS-1 multispectral scanner is capable of providing 
generalized data that could have limited application m urban land 
use studies. The limited spatial resolution and extremely small 
scale of the imagery imposed important limitations on the level 

of information that could be extracted by using conventional 
image-interpretation techniques. The spatial resolution of the 
scanner was not adequate to resolve the details of the many rela- 
tively small objects in an urban scene. 

2. The spectral energies recorded by each picture element 
were integrated measurements the magnitude of which depended on 
the proportion of the picture element occupied by each of the 
objects in the field of view. Consequently, the same spectral 
signature could be derived from a great variety of different 
combinations of surface reflectivities. 

3. In using computer classification programs to classify 
these heterogeneous scenes, serious difficulties were encountered 
in finding spectrally homogeneous urban features of sufficient 
size to be used as training fields. Although clustering tech- 
niques were used to group these heterogeneous picture elements 
into great numbers of similar clusters, the ground-truth meaning 
of these clusters then had to be determined and grouped manually 
into meaningful spatial patterns corresponding to known urban 
land use categories. 

4. The greatest source of classification error was the com- 
puter classification of urban features in which vegetation was a 
major component of the urban scene. For this reason, greater 
classification accuracies could be achieved by making comparative 
analyses of data obtained during different vegetative seasons. 
This statement should at least be true for classification of 
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residential areas, which are normally the most extensive of the 
land use categories within an urban area. 

5. Although some Level II urban categories can be classified 
to varying degrees of accuracy, the current ERTS-1 system and the 
analytical procedures used in this investigation would find the 
most immediate utility for urban and regional planners in pro- 
viding frequent boundary revisions of the urban fringe. The 
greatest spectral contrast occurs between areas of dense vegeta- 
tion and new urban developments. These gross changes in land- 
scape appear most pronounced where forested lands are being 
cleared for urban development. 
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8.0 THE ERTS-1 SIGNATURE EXTENSION ANALYSIS 


8.1 OBJECTIVES 

The objective of the signature extension analysis was to 
measure the spectral signature of a ground feature and use the 
measured signature to identify a similar ground feature in another 
location or at another time. The simplest form of such a proce- 
dure was to use constant signatures that depended only on ground 
features. One objective of the signature extension analysis was 
to test the usefulness of such constant signatures. 

The second major objective was to identify the sources of 
variability in the ERTS-1 data to understand how the signatures 
changed from scene to scene for a constant, unchanging target. 
Freshwater was selected as the test feature because of its homo- 
geneity over large areas and its invariability over periods of 
time. Five water bodies were selected for ground-truth data 
acquisition, statistical training fields, and test sites. These 
bodies were Lake Houston, Lake Livingston, Sheldon Reservoir, 

Lake Somerville, and Steinhagen Lake (fig. 8-1). 

8.2 ANALYTICAL APPROACH 

The feasibility of extending spatial and feature classifica- 
tion was evaluated by analyzing 1972 ERTS-1 MSS data of the 
Houston Area Test Site. Atmospheric haze and Sun angle were 
evaluated because they can affect feature classification. The 
standard data set was the ERTS-1 MSS data for August 29, 1972, 
of Lake Houston, Lake Livingston, and Sheldon Reservoir. Exten- 
sion data sets included the standard data set plus August 28 and 
29 and October 3 and 4 data for Steinhagen Lake; October a ard 
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November 27 data for Lake Livingston, Lake Houston, and Sheldon 
Reservoir; and July 25 and August 30 data for Lake Somerville. 

Several degrees of extension in both time and distance were 
attempted, and the reasons for changes in classification perform- 
ance were analyzed. Extensions from 1 day to 90 days were tried, 
and distances of nearly 160 kilometers {100 statute miles) could 
be evaluated within one ERTS-1 frame. 

8.3 RESULTS 

The following results were obtained. 

1 . The largest feature-dependent variable affecting spatial 
and short-term extension was water turbidity measured in parts per 
million of suspended solids. This parameter varied from 2 to 

5 p/m in Lake Livingston and Lake Somerville to 90 p/m in Lake 
Houston and Steinhagen Lake. 

2. Application of a semiparametric , untrained, discriminant 
technique (ISOCLS) to the ERTS-1 MSS data resulted in seven 
classes of water. These classes described the areas of Lakes 
Livingston and Houston and Sheldon Reservoir more than 0.6 meter 

(2 feet) deep. 

3. Figure 8-2 is an example of the output obtained from 
ISOCLS for Lake Houston and Sheldon Reservoir. The shades of 
gray in the main portions of both sites indicate the five classi- 
fications of water that were obtained from these sites as a func- 
tion of turbidity. Two more classes of water were similarly 
obtained from Lake Livingston. (Statistics are provided in 
table 8-1 . ) 

4. Seven other classifications of water also were obtained 
from these three lakes as functions of shallow water, vegetation 
in the water, and the ratio of water to land in the picture 
element. 



Figure 8-2.- The ISOCLS output 
for Lake Houston and Sheldon 
Reservoir. 
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5. The signatures of the five lakes were extended for the 
same day. The same lake signatures were extended for 36 and 
90 days using a maximum-likelihood (LARSYS) technique. This is a 
parametric, trained classification method that includes the 
statistical means, variances, and covariances that describe each 
class. The normal operation of LARSYS involves selecting train- 
ing fields from which the statistics are calculated. This 
approach did not work well, however, because training fields 
could not be found that yielded unimodal statistics and at the 
same time represented the various types of water in the lakes 
(fig. 8-3). Therefore, the types of water were determined using 
the previously mentioned ISOCLS technique and these statistics 
then were used in LARSYS to classify the water of the other lakes 
and subsequent data sets (fig. 8-4) . With one exception, use of 
this type of extension (ranging from same-day coverage by ERTS-1 

TABLE 8-1.- STATISTICS DERIVED FOR SEVEN TYPES 
OF WATER FROM ISOCLS 
[MSS units, Aug. 29, 1972] 


Water 

type 

Mean value of 
MSS channel no. 


Turbidity 

level 

Site 

1 

2 

3 

m 



1 

35.3 

29.7 

17.1 

3.2 

Highest 

Lake Houston (West 
Fork of San 
Jacinto River) 

2 

33.1 

26.9 

13.9 

2.5 

— 

Lake Houston 

3 

32.1 

24.7 

14.4 

2.5 

— 

Lake Houston 

4 

31.8 

24.4 

12.6 

2.4 

— 

Lake Houston 

5 

26.6 

18.2 

12.9 

2.9 

— 

Sheldon Reservoir 
and Lake Houston 

6 

23.5 

12.9 


1.8 

— 

Lake Livingston 

7 

23.2 

12.8 


m 

Lowe s t 

Lake Livingston 
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to 36 days later) revealed that the variations in atmospheric 
haze were insignificant in water classification, especially when 
considered with respect to the changes in the water caused by 
recent rain and wind disturbance. The exception was a relatively 
thick cirrus cloud that covered the western portions of Lake 
Somerville August 30, 1972. This cloud raised the radiance levels 
of that portion of the lake as much as 10 MSS units. Another 
possible exception, which was unverif iable , occurred on the 
August 28 and 29 coverages of Steinhagen Lake. The radiance of 
the lake was changed approximately three ERTS-1 MSS units over a 
1-day period (table 8-IX) . At that time, no ground truth was 
being acquired from this site. Therefore, there was no way to be 
sure whether this change was due to atmosphere haze or some 
physical change such as increased wind. A ground-truth effort 
was established at this site, but the phenomenon was not recorded 
again. Special-purpose computer programs were developed to aug- 
ment the available software in performing specific signature 
extension investigations. 

TABLE 8-II . - STATISTICAL MEANS DERIVED USING ISOCLS 
FOR STEINHAGEN LAKE 3 


Water type 

Aug . 28, 1972, 
channel no. 

Aug. 29, 1972, 
channel no. 


1 

2 

3 

n 

1 

2 

3 

n 

Least turbid 

23.4 


9.0 

1.5 

26.0 

17.9 

11.3 

2.3 

Most turbid 

34.5 




13.2 

2.2 

37.0 

30.8 

16.6 

3.3 


a Possible haze effect at same locations in lake. 
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6. The MSS data tapes received from G3FC had a cyclic 
striping related to the MSS sensors that was very evident when 
the data were displayed. A computer program was developed that 
eliminated this condition without affecting the data statistics. 

7. A computer program was developed that made it possible 
to merge two strips of MSS data edge to edge, a significant accom- 
plishment with regard to signature extension because in the 
satellite groundtrack, Lake Livingston always appears in two 
computer-compatible tapes. 

8. A computer program was developed that enabled selection 
of specified picture elements and output of their respective data 
values in the four channels of the MSS. 

9. A technique was developed that is much faster than 
LARSYS and ISOCLS for finding water in an ERTS image. In the 
linear discriminant technique for locating water, the data levels 
of only the first and last channels of a picture element from the 
MSS data are evaluated and a determination of whether that element 
is water is made. This technique was later used by the State of 
Texas for their water impoundment study for the National Program 
of Inspection of Dams. 

10. The ISOCLS was modified to include a larger number of 
classes in the low-radiance region. This modification was neces- 
sary for a better evaluation of water, which has relatively low 
reflectances in all four channels. 

8.4 CONCLUSIONS 

The following conclusions can be stated. 

1 . A capability to do short-term temporal (same day to 
36 days) and moderately long-term spatial (within and between 
three ERTS-1 MSS frames) signature extension has been verified 
with respect to large, relatively homogeneous features. 
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2. Ninety-day temporal signature extension for a single 
lake was degraded by the change of Sun angle. The lower Sun 
angle of late fall and winter caused the radiance levels of the 
five sites to decrease by as much as 10 MSS units, even in the 
green band, in which random changes are usually 1 to 2 units 
(table 8-III). No attempt was made to compensate for this type 
of change. Long-term signature extension would require modifica- 
tion of the ERTS-1 MSS data for significant changes of Sun angle. 
Therefore, a data bank of spectral signatures would have to be 
developed on a seasonal basis. 

3. Normally occurring variations in atmospheric haze con- 
ditions appeared to have no major effect on the feature signatures 
in this study. 

4. Haze changes the absolute signature significantly; how- 
ever, in this investigation, it was always by the same amount and 
thus the measured signature of water was always the same. 

TABLE 8-III . - THE MSS CHANNEL 1 (GREEN BAND) DATA LEVELS 

IN MSS UNITS 


Site 

Aug . 29, 1972 

Feb. 25, 1973 

Lake Houston (West Fork 
of San Jacinto River) 

35.3 

27.0 

South-central Lake 
Houston 

33.1 

28.0 

Inlet, northeast 

Lake Houston (cleanest 
water in lake) 

26.6 

21.0 
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5. The greatest difficulty in extending the signature of 
freshwater is caused by turbidity. To recognize a turbid body 
of water, the signature data must be derived from another body of 
water with similar turbidity. 

Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, November 20, 1974 
641-14-07-50-72 
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APPENDIX B 

GLOSSARY OF TERMS, ABBREVIATIONS, 
AND ACRONYMS 
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across- track — across the direction of the spacecraft groundtrack, 
sometimes called horizontal when referring to output product 
coordinates 

ACV - additive color viewer, a device enabling the color enhance- 
ment of one or more black-and-white images of the same scene 
by film density slicing and/or additive color procedures 

ACVP - additive color viewer printer, an ACV that incorporates a 
photographic printing device 

ADP - automatic data processing, such as ccmputer-aided 
computation 

along-track - in the direction of the spacecraft groundtrack, 

sometimes called vertical when referring to output product 
coordinates 

ASCS • Agricultural Stabilization and Conservation Service, an 
agency of the U.S. Department of Agriculture 

CCT - computer-compatible tape containing digital ERTS-1 data. 

The tapes are standard 1.27-centimeter {0.5 inch) wide mag- 
netic tapes in 9-track or 7-track format. Four CCT's are 
required for the 4-band multispectral digital data corre- 
sponding to one scene in the ERTS-1 images. 

clustering - a mathematical procedure for organizing multispec- 
tral data into spectrally homogeneous groups. Clusters 
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require identif ication and in terpretation in a postprocess— 
ing analysis. Both ISOCLS and NSCLAS are spectral cluster- 
ing programs. 

contrast - the ratio of two adjacent scene radiances expressed as 
a number equal to or greater than one 

DAS - data analysis station, a computer system consisting of tape 
drives and computer, a display and control console, and a 
film recorder. The DAS is used to reformat, analyze, and 
review digital remotely sensed data. 

EOD - Earth Observations Division of the NASA Lyndon B. Johnson 
Space Center, Houston, Texas 

EREP - Earth resources experiment package consisting of the 
Earth resources remote sensors mounted on the Skylab 
spacecraft 

ERIPS - Earth Resources Interactive Processing System, a JSC 

system that allows real-time interaction by an investigator 
with several digital spectral analysis procedures. Major 
subsystems include pattern recognition by maximum-likelihood 
classification, image registration, image composition, image 
manipulation, and display. 

ERTS-1 - the first Earth Resources Technology Satellite. The 

ERTS-1 was launched into a circular. Sun-synchronous, near- 
polar orbit at an altitude of approximately 915 kilometers 
(494 nautical miles) in June 1972. It orbits the Earth 14 
times a day and views the same scene every 18 days. 
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geometric accuracy, geographic (latitude-longitude) - based on 

the standard Earth-fixed coordinate reference system, which 
employs latitude and longitude 

geometric accuracy, positional - the capability to locate a point 
in an image with respect to a map 

grayscale - a scale of gray tones between white and black broken 
into an arbitrary number of segments. The ERTS-1 images 
have a 15-step grayscale exposed on every frame of imagery. 
The scale gives the relationship between the gray level on 
the image and the electron beam density used to expose the 
original image. 

ground-control point - any point that has a known location on the 
Earth surface and can be identified in ERTS imagery 

GSFC - NASA Goddard Space Flight Center, Greenbelt, Maryland 

GSFC color composite - color composite of three channels of 

ERTS-1 MSS digital data supplied to users by GSFC. They 
are third- or fourth-generation images as compared with 
first-generation composites produced from CCT's using a 
film recorder. 

HATS - Houston Area Test Site encompassing 18 counties in south- 
east Texas. Houston is the primary urban area in the test 
site . 

image skew - image distortion caused when the scan of the sensor 
is not perpendicular to the plane formed by the spacecraft 
and the instantaneous groundtrack velocity vector 

IR - infrared 


irradiance - the amount of energy impinging upon a unit normal 

surface, per unit time, per unit wavelength, per unit solid 
angle 

ISOCLS - Iterative Self-Organizing Clustering System, a computer 
program developed at JSC that uses a clustering algorithm to 
group homogeneous spectral data. Several controlling inputs 
allow investigators to control the size and number of clus- 
ters. Because the system produces a classification-type 
clustering map in which clusters require postprocessing 
identification and interpretation, the system is frequently 
referred to as a nonsupervised classification system. 

JSC - NASA Lyndon B. Johnson Space Center, Houston, Texas 

A - lambda, the Greek symbol used to designate a wavelength of 
the electromagnetic spectrum 

LARSYS - the name designating the set of classification programs 
developed at the Laboratory for the Applications of Remote 
Sensing (LARS) at Purdue University 

maximum-likelihood ratio - the maximum-likelihood ratio in remote 

sensing is a probability decision rule used to classify a 

target from multispectral data. Two types of errors are 

feasible: failure to classify the target correctly and 

misclassif ication of background as the target. In its 

simplest form, the likelihood ratio is PVP, or the prob- 

t: d 

ability P of an unknown spectral measurement being 
classified as target t, to the probability of an unknown 
spectral measurement being classified as background b. 

When p t / p j-, 1 1 > the decision is t; when p t / p j-, < 1 > the 
decision is b. Probability density functions are computed 
from spectral samples, often referred to as training 
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samples. As the number of training samples increases, the 
mathematical computations of the maximum- likelihood ratio 
increase in complexity. As a result, digital computer 
analysis is required, and the entire process is referred to 
as automatic data processing (ADP) of multispectral remotely 
sensed data, or automatic spectral pattern recognition of 
multispectral remotely sensed data. 

MCFV - multichannel film viewer, a device used to electro- 

optically enhance several bands of black-and-white imagery of 
the same scene. Enhancement is achieved through film density 
slicing and/or additive color viewing. 

MSDS - multispectral data system, which includes an aircraft 

24-channel scanner and a ground data analysis station. The 
latter is one of the two major data analysis stations in the 
EOD/JSC Data Analysis Laboratory. (See DAS.) 

MSS - multispectral scanner system, sometimes referred to simply 
as the multispectral scanner; usually refers to the opera- 
tional scanning system on ERTS-1 

MTFO - module training-field option, a supervised computer-aided 
technique using the LARSYS classifier to allow modification 
of the statistical inputs 

multispectral scanner spectral bands - the division of the visible 
and near-infrared portions of the electromagnetic spectrum 
into discrete segments 
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nautical mile - equivalent to 1/60° at the Earth Equator, or 
approximately 1852 meters (6076 feet) 

nonsupervised classification - a procedure by which spectral data 
are grouped into homogeneous clusters. Identification and 
interpretation are achieved in a postprocessing analysis. 

NSCLAS - the name of a clustering computer program developed by 
the Laboratory for the Application of Remote Sensing at 
Purdue University (See clustering.) 

pixel - picture element, refers to one instantaneous field of 

view (IFOV) as recorded by the mul tispectral scanning sys- 
tem. On the ERTS-1 system, a pixel is equivalent to approx- 
imately 0.44 square hectometer (1.09 acres). One ERTS-1 
i frame contains approximately 7.36 * 10 6 pixels, each 

described by four radiance values. 

radiance - a measure of the radiant energy emitted by a radiator 
in a given direction 

reflectance - the ratio of the radiance of the energy reflected 

from a body to that incident upon it; commonly measured in 
j percent 

i 

i 

1 
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scene registration - the capability to superimpose points on two 
images of a scene taken at the same time 

SHNF - Sam Houston National Forest, located in the Houston Area 
Test Site 

spectral response - the spectral radiance of an object sensed at 
the satellite and recorded by the multispectral scanner 
system 


Sun azimuth angle - angle in degrees measured in the horizon 

plane from true north to a vertical circle passing through 
the Sun 

Sun elevation angle - angle of the Sun above the horizon measured 
in degrees 

supervised classification - a classification procedure in which 
data of known classes are used to establish the decision 
logic from which unknown data are assigned to the classes. 

The ADP supervised classification procedure used at JSC 
during the ERTS-1 project incorporated a Gaussian maximum- 
likelihood decision rule. 

swath path - the dimension on the ground scene transverse to 
spacecraft velocity and within the sensor field of view 

system-corrected images - film images generated by a data- 

processing subsystem that makes initial radiometric and geo- 
metric corrections as the video-to-f ilm conversion images are 
recorded on 70-millimeter film through an electron beam 
recorder; formerly referred to as bulk images 
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temporal - that which exists in the physical world as it relates 
to time 

temporal registration - the capability to superimpose two images 

of the same scene taken at different times (same or different 
spectral bands) 

test field - the spatial sample of digital data of a known ground 
feature selected by the investigator used to validate the 
statistical parameters generated from training field samples 

threshold - the boundary in spectral space beyond which a data 
point (pixel) has a sufficiently low probability of being 
included in a given class and, therefore, is purposely 
excluded from that class 

training field - the spatial sample of digital data of a known 
ground feature selected by the investigator from which the 
spectral characteristics are computed for use in supervised 
multispectral classification of remotely sensed data. The 
statistics associated with training fields form the input 
to the maximum likelihood ratio computations and, in a sense, 
"train" the computer to discriminate between samples. 

USDA - U.S. Depax cment of Agriculture 

USGS - U.S. Geological Survey of the Department of the Interior 

UTM grid - Universal Transverse Mercator grid, a rectangular coor- 
dinate system derived from a transverse mercator projection 
by which points or areas on the Earth surface car. be readily 
described and located within a unique quadrilateral area on a 
map. Precision processing of the ERTS-1 imagery corrects the 
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imagery scene by scene so that its geographic orientation 
will conform to this UTM grid system. 
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